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ABSTRACT 


A fluctuating extended radar target is a general model for 
all classes of clutter-producing media. The complex reflectivity 
of a target can be estimated. Based on this data, and under 
the assumpti on of WSSUS, various properties of interest can be 
investigated. Using common tools of analysis, three programs, 
have been developed for evaluating: 

1. The probability distribution of clutter samples 
obtained from a single range cell. 

2. Correlation between the clutter data of two 
range cells. 

3. Clutter Power Spectral Density of a single range cell. 

These programs (in FORTRAN IV) were run on mini-computer PDP-11/24 
with generated test data and produced satisfactory results. Some 
limitations of the methods used are discussed and a few 
suggestions for improvement given. 
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CHAPTER 1 


INTRODUCTION 


1.1 PREAMBLE 

When a beam of electromagnetic energy is transmitted from 
a directional radar antenna, the signal received by the radar 
would comprise of 

a) reflection from the desired target (s), if any, 

b) unwanted reflection from other than the desired 
targets in the path of the beam, termed CLUTTER, 

c) external noise generated by electrical disturbances, 

d) thermal noise at the receiver input. 

The difference between clutter and noise is that an increase in 
the transmitted power results in an increase in the clutter power 
as well as the desired signal level but not in the noise level. 

The classical problem in radar signal processing is to 
detect the desired target signal which is embedded in clutter. 

The obvious method of doing this is to enhance the target signal 
relative to the clutter by exploiting the differences in spatial, 
temporal and spectral properties of the target as compared to the 
clutter source. 

For example the traditional method of clutter suppression 
is the MTI (Moving Target Indication) filter. This works on the 



assumption that the target of interest is moving radially to 

source 

the radar while the clutter/is not. This assumption is not 
always correct. 

More recently, various adaptive clutter suppression 
filters have been designed or theoretically formulated. Their 
principle is to adapt to the prevalent clutter conditions and 
in subsequent processing enhance the signal-to-clutter ratio, 
one way being to subtract an estimate of the clutter out of 
the total received signal, the desired target signal remaining 
basically unaffected. 

A study of the statistical properties of land, sea and 
wea+ clutter is of interest. Extensive studies of clutter 
collection and analysis of data have already been carried 
out in various countries. However, the studies on land clutter 
do not relate to the Indian subcontinent and are, therefore, of 
no particular use to us [l]. 

It would, therefore, be gainful to measure clutter data 
from various types of terrain in this country, under different 
climatic conditions, and to carry out studies on the data. Due 
to the random nature of the clutter and the large masses of the 
data involved, the study would have to be based on statistical 
evaluation. 

The time varying nature of clutter, its degree of spatial 
correlation and spectral content are important aspects of 
such a statistical evaluation. Some methods and programs for 



carrying out this evaluation are considered in this thesis. 
However, these, at most form a preliminary investigation into 
the clutter data and more detailed work in this direction is 
recommended. 

A better understanding of clutter properties would in 
future influence the design of radar receivers/adaptive clutter 
suppression systems. 

1.2 RADAR CLUTTER DATA MEASUREMENT 

1.2.1 Brief Background: 

A clutter signal consists of unwanted radar returns. In 
experimental studies of clutter producing targets, however, the 
clutter is a source of information. Such a study can be viewed as 
consisting of two interlinked parts: 

a) The measurement part 

b) The statistical analysis part 

A vehicle mounted clutter measurement system based on a CW radar 
has been designed and setup at IIT [l]. The specifications are 
given at Table 1.1. It would constitute the source of measured 
clutter data. 

In order to analyse this statistical data, off-line, the 
implementation of relevant statistical methods on a minicomputer 
(PDP-ll/24) using FORTRAN language, was envisaged. This forms the 
subject of the present work. 



The result was the development of 3 programs for stati- 
stical processing of clutter data. These relate, respectively, 
to 

a) Probability distribution 

b) Spatial correlation 

c) Power spectral density 

A consolidated list of environmental and system parameters to be 
recorded 'on-site’ is suggested in Table 1.2, 

1.2,2 Model: 

A general model for a clutter source; is a randomly fluctu- 
ating, range-spread target. Under the assumption of WSSUS (wide - 
sense stationarity and uncorrelated scattering), and using complex 
envelope notation as applied to narrow-band bandpass signals, 
the complex envelope of the received (clutter) signal is given 
by [2] 

s(t) = f VE^ f(t-A ) b( t - X)dX (1.1) 

— KX) 

where, 

is the energy of the transmitted waveform 

#s» 

f(t) is the complex envelope of the transmitted signal 
normalised to unit energy 
X is the two-way range-delay in time units 
b(t,X)is the complex reflectivity of a dX wide strip located 

at range-delay X , 





It should be noted that b includes path losses, antenna gains, 
effect of polarisation, etc* 

The clutter source is equivalent to a random, linear, 
time-varying filter with input VE^ f(t) (complex envelope of the 
transmitted signal) , impulse response b(t - X ) = b(t, X), and 
output s(t) (complex envelope of the received clutter signal). 

The measurement system is designed to provide estimates of 
the samples of the complex reflectivity (or impulse response) b 
as a function of both discrete range and discrete time. 

Integer variables K and n will be used to denote the 
discrete range and the discrete time respectively* 

For a particular set of parameters in Table 1 * 2 ,- the data 
would consist of the equispaced-in-time sequence of sample esti- 
mates bj, * Each bj^ consists of the in-phase and quadrature 
components, i.e* two real numbers* 

If the beam direction is kept fixed and assuming that the 
illuminated patch consists of r range cells given by 


K = K 3 , ..., Kj, ; 0 < K < 32767 (1.2) 

the recorded data would consist of N (e.g* 1024) windows, the 
nth window comprising the sequence 


tfnrnrfn 

^Ki » “ 1,2,...,N ; 


b./‘ = (b^"'^ , bq'^j, ) 






(1.3) 


where I and Q denote in-phase and quadrature-phase respectively* 
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It may be noted that the window index n is the same as 

the discrete time index. for data samples pertaining to a given range 
cell Ki. 

The data will be recorded on a storage device (Mag tape). 

1.2.3 Preprocessing: 

A re-arrangement of the above database is required prior 
to subjecting it to the envisaged statistical processing. 

The original database can be thought of as an Nxr matrix 

(e.g. N = 512, r = 173) whose N rows are the N windows. Let the 

matrix be denoted by [b„ ^]. 

n j K 

The data is now required range-cell-wise. This can be 
achieved by a transposition of the above time-range matrix, 
yielding 

This preprocessing manipulation has been presumed. 

However, it can be included as the first step in the statistical 
processing. . 

The methods in Chapters 2&3 process N-point data sequences 
pertaining to individual cells, i.e. rows of the transposed matrix. 

1.3 AN OUTLINE OF THE THESIS 

The statistical processing/analysis scheme comprises three 
parts, one for each of the properties to be investigated. 

To observe whether the time-variations of the clutter in a 
given range cell match a known distribution, Chapter 2 describes 



qualitative tests for Rayleigh, Vi/eibull and Lognormal models. 

Two examples j llustrate the use of this method. 

The first part of Chapter 3 deals with the estimation of 
clutter power spectral density. A scattering function description 
of a doubly-spread .radar target, based on the assumption of 
WSSUS, is given. An approximate model with a tapped-delay-line 
structure is described. The FFT algorithm used in the program 
is explained. One example is given. 

The second part of Chapter 3 deals with the estimation of 
correlation between the clutter returns from two range cells. 

Two methods (Direct and FFT) for the estimation are outlined and 
some illustrative examples given. 

Chapter 4 summarises the programs developed and gives 
few suggestions for improvement. 

All examples are based on generated data. 



TABLE 1.1 
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MEASUREMENT SYSTEM . SPECIFICATIONS 

System : Truck-mounted S and X Band ON radar 

Transmitted Power : 400W 

IF : 70 MHz 

Modulating Signal ' : PRBS generated by 15-stage SR 

Type of modulation; Biphase 

Clock rates , ; 1.25, 2.5, 5, 10 MHz 

Antenna ; Parabolic dish, dia 1 m at S-band 

m at X-band 

Polarization ; Horizontal/Vertical (switchable) 

Sidelobe level ; -20 dB max. 

Antenna gain ; 25 dB at S-band 

33 dB at X-band 

Scan rate ; 10 rpm, in Azimuth, over 360° 

Angle of depression; . 0.5° to 17.5° (below horizontal) 

Max. unambiguous range corresponds to the 32767th range cell 
Positional accuracy; Azimuth + 0.5° 

Elevation (depression) + 0.25° 

System output ; Samples of complex reflectivity (in- and 

quadrature-phase components) of 15 adjacent 
range-cells are simultaneously estimated, 

(The measurement system was designed at A.C.E.S,, IITK) 



TABLE 1.2 


ENVIRONMENTAL AND SYSTEM PARAMETERS 



Parameter 

Typical example 

1. 

Date 

150688 

2. 

Place 

KANPUR 

3. 

Terrain 

CITY 

4. 

Weather 

CLEAR 

5. 

Temperature (°C) 

28.5 

6. 

Pressure (mm) 

180 

7. 

Humidity {"/*) 

80.5 

8. 

Wind Velocity (km/h) 

0 

9. 

Carrier frequency (GHz) 

3.0 

10. 

Clock-rate (MHz) 

10.0 

11. 

Depression angle (®) 

15.5 

12. 

Azimuthal angle (°) 

250 

13. 

Polarisation (v/H) 

H 

14. 

Time 

1615 

15. 

Samples/ cell 

1024 


CHAPTER 2 


STATISTICAL DISTRIBUTION OF THE CLUTTER SAMPLES FROM ONE 
RANGE CELL: A QUALITATIVE EVALUATION 


2.1 PURPOSE 

For any class of clutter sources, a general model is a 
randomly fluctuating, range-extended (doubly spread) radar target, 
e.g. a thickly vegetated hillside in windy weather. 

The prevailing weather contributes its own factors (such 
as rain, dust, snow, hail, wind) to the random nature of the 
received clutter signal. 

For instance, the spectral content of clutter from a wooded 
area would vary with the wind speed through that area. 

It is of interest to study the clutter returns from various 
types of reflecting surfaces such as cultivated land, water 
stretches, deserts, forests, hills, snow-clad mountains, built-up 
areas, and also to relate the nature of the clutter variations 
to the general class of the source, its physical spread, as well 
as geographical and meteorological factors. 

The extent or spread of the target' (in the illuminated 
space) is obtained by the number of range cells which contribute 
to the clutter returns. 

Estimates of the complex reflectivity of each range cell 
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in the illuminated patch as a function of (discrete) time, 
n 

bj, , where 

b denotes complex reflectivity, 
n denotes the discrete time index, 1 £ n ;< N, 

K denotes the discrete range index. 


can be obtained with the radar system described in Section 1.2» 


For the Kth range cell, if a set of N observations are 
taken, the resulting data comprises of a set of N complex numbers 


u 1 u 2 ? 3 

bj^ , bj, # bj, , 



representing N complex-reflectivity estimates equi-spaced-in-time. 

The question that logically follows is that does the 
distribution of these sample values follow any of the known 
probability distributions? 

Three statistical models have been suggested for the 
magnitude of the complex reflectivity, literature 

(e.g. [ 2 ], [3]) for all types of terrain clutter* These are the 
Rayleigh, Weibull and Lognormal distributions. Their densityy^ 
distribution functions are given in Table 2.1. 

The adequacy of these distributions, in relation to the 
Indian terrain-and-meteorological conditions, for modelling the 
time-varying clutter return, requires verification. One approach 
in carrying out such a check is to assume a particular statistical 
model and then to test the validity of the assumption with the 
data. This is the approach used in this work. 



TABLE 2.1 


RALEIGH, 

P(r) 

F(r) 

P(x) 

F(x) 

p(y) 

F(Y) 

Symbols ; 


,WEIBULL AND LOGNORMAL DISTRIBUTION FUNCTIONS 


RAYLEIGH 


0, r < 0 

1 r ^-r^/2a^ 


= < 


7 " 


0, r < 0 


, r > 0 


1 - e 


, r > 0 
WEIBULL 


> -CO < or < CO 


.1 


0, X < 0 

5 exp[-(|)’^], x>0 

0, X < 0 

1 , X > 0 


T)>0, a>0 


LOGNORMAL 


.1 


0, y < 0 


1 . . .. e- . Y > 0 

721 oy 20^ 


\ a > 0 

<H<oo 


0, y < 0 

Y 1 dy, y > cl 

^ ® 20^ i 

p: probability density function (pdf), F: Cumulative 
Distribution Function (CDF)> r,x,y: random variables (RV) 
Y: particular value of y. 


13 


2.2 TEST FOR RAYLEIGH AND WEIBULL DISTRIBUTIONS 

The data relates to the random variable 

X = b = Ibj (2.1) 

for a given range-cell, whose distribution is to be studied. It 
suffices to test for a Weibull distribution only, the Rayleigh 
model being a special case of the Weibull model. The Rayleigh 
model shall be valid when the following two conditions are met: 

a) Weibull model holds good 

b) the Weibull parameter r) = 2 

A simple method to assess the validity of an assumed (in 
this case Weibull) distribution is to manipulate the measured 
statistical data in such a manner that, when plotted graphically, 
it: 

a) plots into a straight line if the statistical model 
is valid , 

b) does not plot into a straight line in all other cases. 


The cumulative Distribution Function of the Weibull distribution 
is 


F(x) = 1 - exp [- (x/o) ■], X >_ 0, Ti > 0, a > 0 


W. 

After algebraic manipulation this can written as 


( 2 . 2 ) 




= T) In X - Tj Ino, X > 0, t) > 0, a > O 


(2.3) 
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Fig. 2.1; SCHEME FOR TESTING THE VALIDITY OF RAYLEIGH AND 
WEIBULL MODELS 
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which is the equation of a straight line with slope t) and y-axis 
intercept (-t) In ct) . 

The various steps in the complete test procedure are shown 
in the flow chart of Fig. 2.1. This forms the basis for the 
first half of Program No. 1, the second half consisting of a similar 
test for the lognormal distributional model. 


2.3 TEST FOR LOGNORMAL DISTRIBUTION 

The Cumulative Distribution Function for the Lognoimal 
model is, from Table 2.1, 

= I ]db, (2.4) 

b>_0, 0>O, - 

where, 

the random variable b = |b| 

B is a particular value of b 

\i is the mean of Inb 

(P' is the variance of Inb 

With a change of variable x' = Inb, eqn. (2.4) becomes 
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(in , 2 

F(B) = f exp[- ]dz, o>0, (2.6) 

— oo<2<«>> — oo<|i<oo 

which is the error function in classical Detection and Estimation 
theory. 


F(B) = erf ( 



= 


(2.7) 


Since z = , we put Z = iHizH (2.8) 

where Z is a particular value of the RV z, . 


F(B) = erf(Z) = Q(Z), B>0,-~<Z<co (2.9) 

If the lognormal model is valid, a plot of F(B(I)) against 
Q(Z(l)), I = 1,2,...,N would be a straight line through the origin 
with unity slope. 

On the face of it this test (for the lognormal model) seems 
more rigorous than that (in 2.2) for the Rayleigh and Weibull 
models. However, the calculation of Q(z) = Q( ~0 . ^T il;) requires the 
values of p and o, which are not previously known and have to be 
evaluated from the data itself . Thus p and a are themselves RVs , 
which would vary with each set of measured data. The introduction 
of the two RVs has to be offset against the greater ease in 
(visual) assessment from the graphical plot. The value of the 
error function integral 

Q(Z) = erf(z) = / 7 ^ dz (2.10) 

—00 « 




Fig. 2.2: SCHEME FOR TESTING THE VALIDITY OF LOGNOR\UL MODEL 





is evaluated (by a subroutine) using the Gaussian quadrature 
method of Numerical Analysis [5]. 

The complete procedure involved in this test is shown as 
the flowchart in Fig. 2.2. The first 5 steps are identical to 
those in Fig. 2.1, being, in the program, subroutines commonly 
shared. 

2.4 DISCUSSION AND SUMMARY 

1. Successive valuesof the measured parameter 

b(n), n=l,...,N (2.11) 

should be sufficiently separated in time as to be uncorrelated. In 
case this condition is not satisfied during the measurement phase, 
the entire period corresponding to n=l,...,N can be divided into 
intervals equal to the requisite separation (except for the last 
interval v^ich may be of a smaller length) and the samples averaged 
out in each of these intervals. For this reason a linear time 
averaging facility has been provided in the menu-driven program 
and is optional. (The number of consecutive samples to be averaged 
is set by the user). . 

The number of data points after averaging, if any, should 
be sufficient for meaningful post-processing. This entails taking 
an adequate number of data points (N) to start with. 

2. The testing is qualitative and, in each case, provides for 
one of the following decisions (from the graphical plots) : 


IV 


a) the model under test is adequate 

b) the model under test appears to be adequate 

c) the model under test is rejected as inadequate 

d) the number of points are too few to make a decision. 

As the number of data points is increased, the assessment becomes 
correspondingly easier upto a point, whereafter further increase 
would only cluster the plotted points together unless the scale 
is correspondingly expanded. 

However, in the lognormal test, as the mean and variance 
of Injbj, (^, 0 "^), are estimated from the data itself j the accuracy 
of the estimates increases with N, 

For example, the estimates of (p.,cj^) were (.3,1.7) with 
N = 16 when the test was run with data from a simulated lognormal 
model with actual (p,a ) = (0,1). Increase in N gave more accurate 

o 9 

estimates for (p,cr ), For N = 1024, ( 11,0 ) were (.00278, 1.02), 
the plot was a straight line with slope = 1.0107 ' and the Y-axis 

intercept = -4.758x10'*^ 

3. The basic limitation of this model-fit test procedure is 
its qualitative nature. It is subjective and provides a prelimi- 
nary insight into the time-varying behaviour of the data. It is, 
therefore, suggested that a second test, which is quantitative in 
nature, be included in the testing schemes, at a later stage, for 
corroboration-, especially in cases where the decision is of the 
form 2(b) above. 





4* The program which has been written incorporating the method 
outlined in this chapter is menu-driven and essentially consists of 
two parts : 

a) the Rayleigh/Weibull model validity assessment ■ 

b) the Lognormal model validity assessment 

The user is provided the facility to: 

a) carry out one or both of these two tests 

b) repeat one or both tests 
without having to re-run the program. 

. 5. As an illustration, the results of running the program with 

generated test data, for the following two cases, are shown: 

a) Rayleigh distributed data 

b) Lognormal distributed data 

on subsequent pages. (The data was generated with the. help of the 
program, in .Appendix A.4) , 



Example 1: 


We consider a case where we are given a sequence of 
Rayleigh distributed clutter samples which we subject to 

the two tests for (a) the R-ayleigh/Weib«41 (b) the Lognormal 
mod<l respectively. 

Table 2.2 shows N = 32 samples of the complex reflectivity 
pertaining to a range cell (say Gth cell), whose magnitudes follow 
a Rayleigh distribution with parameter o = 1.5. 

No time-averaging was done on the data. 

The data was processed by the Rayleigh/Weib<ull test program 
segment and the resultant data is shown in Table 2.3 and its 
plot in Fig, 2.3, The straight line obtained indicates that the 
data is in general, Weibokll distributed and the slope t) = 2.0 
indicates that it is in particular, Rayleigh distributed. The 
Rayleigh distribution parameter estimate obtained from the plot 
is 

Og = 1.5000007 

which compares favourably with the a value of the generated data. 

The same data was next processed by the lognormal test program 
segment and the resultant data is shown in Table 2.4 and is plottec 
in Fig. 2.4. The best straight line fit is indicated. The plot 
is not a straight line which appears to indicate that a lognormal 
model for the given data is not adequate. However, as pointed 
out in Section 2.4, the lognormal model test method is based on 
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TABLE 2.2 
Cell d a t a 

no. In-ohase part auadrature part 

1 0.27911732E + 01 0 . 852 102r2E-01 

2 0.55312514E+01 0 . 67858332E+00 

3 0.21003S73E+01 0 . 95&28780E+00 

4 0.49009085E+00 0 . 24491093E+01 

5 -0.92636120E + 00 -0 . ie216.825E + 01 

G 0.20585775E+01 0 . G6408443E+00 

7 0.80bS'123GE + 00 0 . 1G638S13E + 01 

8 -0.17642403E+01 -0 . 808G4424 E+00 

9 -0.G3073218E-01 -0 . 1&882457E+01 

10 0.12802414E+01 0 . 121&G1 16E+01 

11 -0.131&0012E+01 0.81901288E+00 

12 -0.37892115E+00 0 . 15727392E+01 

13 0.14483830E+00 -0 . 1 4 1 &2G35E+01 

14 0.39775985E+00 0 . 14313120E+01 

15 -0.49155420E+00 -0 . 1 209G 145E+0 1 

16 0. 13566177E+01 -0 . 1 4005318E4-00 

17 0.11907479E4-01 0 . 72769538E-0 1 

18 0. 12204448E+01 0 . 26475656E+00 

19 0.82700241E+00 0 . 69927490E+00 

20 -0.66707754E+00 0 . 92172498E+00 

21 0.88461250E■^00 0 . 406438 1 4E + 00 

22 0. 12938941E+00 0 . 1 020 1796E+0 1 

23 -0. 16734710E+00 -0 . 84559977E+00 

24 0.33419365EH-00 0 „ 85520470E + 00 

25 -0.57129407E+00 -0 . 47860053E+00 

26 -0.1421412SE+00 -0 . 791 88430E+00 

27 0.57653236E+00 0 . 2233461 OE+00 

28 0.59902543E+00 -0 . 32917646E+00 

29 -0.469g5527E + 00 -0 . 251 42889,E-0 1" 

30 -0.50128603E+00 0 . 221 71623E+00 

31 0.24894603E+00 0 . 97264990E-01 

32 0.31414589E + 00 0 - 21 569446E + 00 



OUTPUT DATA 


TABLE 2.3 

UEIBULL DISTRIBUTION TEST : 


1 

X 

-0. 13194S65E+01 

Y 

-0.34499032E+01 

o 

-0. 964781 17E+00 

-0 .27404928E+01 

3 

-0.75368863E+00 

-0.23183076E+01 

4 

-0.60124397E+00 

-0 .20134182E+01 

tr* 

U 

-0.48081008E+00 

-0. 17725508E + 01 

b 

-0.38051108E+00 

-0.15719523E+01 

n 

i 

-0.2940014.9E + 00 

-0.13989335E+01 

8 

-0.217484312+00 

-0.12458992E+01 

9 

-0. 14850001E+00 

-0.11079304E+01 

10 

-0.85358232E-01 

“'0-981&4696E + 00 

11 

-0 .26S42147E-01 

-0.86461544E+00 

12 

0.27957678E-01 

-0.75501478E+00 

13 

0.7974749SE-01 

-0.65143549E+00 

14 

0.12908906E+00 

-0.55275208E+00 

15 

0 . 17644542E+00 

-0.45803934E+00 

16 

0 .22220883E+00 

-0.366512891+00 

17 

0 .266721S4E+00 

-0.2774a659E+00 

18 

0.31029534E + 0;) 

-0. 19033924E+00 

19 

0 .35322423E+00 

-0. 10448&90E+00 

20 

0 .39573688E+00 

-0.19356817E-01 

21 

0.43828440E+00 

0.6563S542E-01 

O •“) 

0.48103148E+00 

0.15113252E+00 

23 

0„52438730E+00 

0.237S4402E+00 

24 

0.56878209E+00 

0 .326634 2 9E + 00 

-f-;* tr 

0.61476314E+00 

0.41859576E+00 

26 

0.66306633E+00 

0.51520199E+00 

27 

. 0.71475738E+00 

0.&1S58422E+00 

28 

0.77151549E+00 

0.73209941E+00 

29 

0.83630300E+00 

0.86167556E+00 

30 

0.91535592E+00 

0.10197815E+01 

31 

0. 10269278E+01 

0 .12429250E + 01 

32 

0. 17173835E+01 

0.26248367E+01 





TABLE 2.4 

LOGNORMAL DISTRIBUTION TEST ; OUTPUT DATA 


X Y 


1 

0.71114302E-02 

0.31250000E-01 

o 

jut 

0 .30459225E-01 

0.62500000E-01 

3 

0 .&2957972E-01 

0.93750000E-Oi 

4 

0 .998S4063E-01 

0. 12500000E+00 

5 

0.13870990E+00 

0. 15625000E + 00 

6 

0.17804623E+00 

0. 18750000E+00 

7 

0.21711010E+00 

0.21875000E+00 

8 

0.25545776E+00 

0.25000000E+00 

9 

0 .292S4254E + 00 

0.28125000E+00 

10 

0 .329136&7E+00 

0.31250000E+00 

11 

0.36428505E+00 

0.34375000E+00 

12 

0.39827728E+00 

. 0,37500000E+00 

13 

0.43113327E+00 

0.40625000E+00 

14 

0.46288955E+00 

0.43750000E+00 

15 

. 0.49359494E+00 

0 .46875000E+00 

16 

0.52330554E+00 

0.50000000E+00 

17 

0.55208200E+00 

0.53125000E+00 

18 

0.57998836E+00 

0.56250000E+00 

19 

0 . 6 0 7 0 9 5 9 0 E + 0 0 

0.59375000E+00 

20 

0.63346857E+00 

0 .62500000E+00 

21 

0.&591S660E+00 

0.65625000E+00 

22 

0.68432975E+00 

0.68750000E+00 

23 

0.70898724E+00 

0.71875000E+G0 

24 

0.73326164E+00 

0.75000000E+00 

25 

0.75727671E+00 

0.78125000E+00 

26 

0.78118807E+00 

0.81250000E+00 

27 

0.80520743E+00 

0.84375000E+00 

28 

0.82964867E+00 

0-S7500000E+00 

29 

0.85503119E+00 

0.90625000E+00 

30 

0.88237929E+00 

0.93750000E+00 

31 

0.91444021E+00 

0.96875000E+00 

32 

0 .99367672E+00 

0. 99999899E+00 
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the estimation of the mean and variance of the In of the data 
magnitudes from the data itself. Therefore, a large number of 
sample values (e.g. N = 1024) should be taken in order to obtain 
a correct assessment (when subjecting the cell data to the 
lognormal test program segment). This is further illustrated in 
the following example. 

Example 2: 

Here we consider a sequence of N = 1024 samples of the 
complex reflectivity of a particular range cellH. The 
magnitude values of the complex samples follow a lognormal 
distribution (within an absolute error of + 1.0 x 10 ^ due to 
the data generation method). The distributional parameters 
are ; 

mean of In [bj^j = p = 0.0 

I I 2 

variance of In |bp^j = cr = 1.0 

Table 2.5 shows 32 of these 1024 sample values beginning with 
the 32nd sample followed by samples at equispaced intervals 
of 32 samples. 

No. time-averaging was done on the data. 

The data was first assumed to satisfy a Weib^ull distribution. 
The 1024 point data was processed by the Rayleigh/WeibUll test 
program segment. The resultant data is shown* in Table 2.6 and 
its plot in Fig. 2.5. We conclude from the plot that the 
Weibull model does not fit the data. 





Next the data was assumed to satisfy a lognormal distri- 
bution, The 1024 point data was processed by the lognormal test 
program segment. The resultant data is shown* in Table 2.7 and 
its plot in Fig. 2.6. The plot is a straight line with 

slope = 1.0107476 

Y-axis intercept = -0.475845x10"^ 

indicating the model to be an adequate description of the 

2 

clutter data. The estimates of p and o obtained from the data 
are: 

Pjc = 0.00278 

= 1.02 

These may be compared with the corresponding estimates obtained 
when using only the N = 32 data values of Table 2,5: 

Pg = a. 1417 

= 1.4145 
h 


only 32 output values 


out of 1024, beginning with the 32nd 


value, followed at intervals of 32 values. 



T A El L b b . D 
CELL DATA 


IN~FHASE PAET 

1 0.643S2300E+01 

2 0.9244&465E+02 

3 0.340014&7E+01 

4- 0.90931777E+00 

5 -0. 1244509 lE+01 

6 0.30067203E+01 

7 0.94849020E+00 

8 -0.22073731E+01 

9 -0.&6622227E-01 

10 0.14229839E+01 

11 -0.12694244E+01 

12 -0.38188795E+00 

13 0. 128972LlE-f-00 

14 0 .3S822775E+00 

15 -0.40718207E+00 

16 0.11641711E+01 

17 0.922S6336E+00 

18 0.97726899E+00 

19 0.b0236237E+00 

20 -0.5009504faE+00 

21 0.60773802E+00 

22 0.91490082E-01 

23 -0. 10879235E+00 

24 0.22325191E+00 

25 -0.35265395E+00 

26 -0.90000726E-01 

27 0.33962879E+00 

28 0.36092383E+00 

29 -0.2&728570E+00 

30 -0.28947607E+00 

31 0. 14460196E+00 

32 0, 17777355E+00 


QUADRATUKE PARI 
0 . 19654922E+00 
0. 11341490E+02 
0. 15480?14E+0i 
0.45470920E+01 
-0.24473181E+01 
0.9S994948E+00 
0.19558196E+01 
-0. 10117555E+01 
-0 . 17832400E+0i 
0.13522596E+01 
0 .7g002595E+00 
0.15S49701E+01 
-0.126ir200E + 01 
0 . 13250426E+01 
-0 .. 1Q019 920E + 01 
-0.12018557E400 
0.563g8459E-01 
0.21200335E+00 
0.50933009E+00 
0.6921 8 123E+00 
0.27922726E+00 
0.72135979E+00 
-0 . 54972440E+00 
0.57130373E+00 
-0.29543519E+00 
-0 , 50140363E+00 
0. 13157070E+00 
■' 0 . 1 9 8 3 3 4 8 7 E + 0 0 
-0. 14299947E-01 
0. 12a03377E + 00 
0.564970i9E-01 
0 . 1 2 2 0 6 0 3 8 E + 0 0 



c-j ITJ kD tv CO CTt O C’‘3 m in 


lEIBULL DISTRIBUTION TEST 


OUTPUT DATA 


X 


Y 


7 

3 

3 

) 

L 

■> 

3 

5 

> 


-0. i8&27377E+01 
-0. 15341281E+01 
-0. 13180081E+01 
-0. 115034S2E+01 
-0. i0099384E+0i 
-0.S8714826E+00 
-0.77641833E+00 
-0.67448S19E+00 
-0.57912910E+00 
-0.4887780SE+00 
-0.40224794E+00 
-0.31863827E+00 
-0.23720007E+00 
-0. 15731005E+00 
-0.7841017SE-01 
0 .OOOOOOOOE+00 
0 .7841008gE-01 
0. i5731007E+00 
0 .23720244E+00 
0 . 31S64005E+00 
0 .40225020E+00 
0 > 4 8 8 7 7 8 2 3 E + 0 0 
0 .57912934E+00 
0 »67448950E + 00 
0 „77641940E+00 
0 .887i4939E+00 


-0.34499032E+01 
-0 .27404928E+01 
-0 .2318307&E+01 
-0 ,20134182E-5-01 
-0 .17725508E+01 
-0 .15719523E+01 
~0.139B9335E+01 
-0.1245S992E+01 
-0.11079304E+01 
-0.981S4696E+00 
-0.864&1544E+00 
-0.75501478E+00 
-Q„&5?.43549E+00 
-0.55275208E+00 
-0.45803934E-f-00 
-0.3fa651289E+00 
-0.27748659E+00 
-0. 19033g24E+00 
-0. 10448690E+00 
-0. I935i3817E-01 
0.65&38542E-01 
0, 15113252E+00 
0.23784402E+00 
0.32663429E+00 
0,4185g576E+00 
0.51520199E+00 


0.10099S9&E-^01 0.61358422E + 00 

0 ,11503501E-!-0i 0.73209941E + 00 


0.13180094E+01 
0 , 15341 194E-f01 
0 . 1S627195E+01 
0. 45340991 E+01 


0.8&167556E+00 
0. 10197815E+01 
0. 12429250E+01 
0.2b24a3&7Er01 
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TABLE 

LOGNORMAL DISTRIBUTION TEST : OUTPUT DATA 


X 


Y 


1 0.33704311E-01 

2 0 .6533444S'E-01 

3 0.S7678g5SE-01 

4 0.i2912902E+00 

5 0,lS037679E + .00 

6 0.19147357E+00 

7 0.22245741E+00 

8 0.25334883E+00 

9 0.2841&g5eE+00 
LO 0 .31493151E + 00 
11 0 .345&5195E+00 
.2 0.37S33f378E + 00 
.3 0.40699780E+00 
.4 0 .437&4171E + 00 
.5 0.46S27841E+00 
S 0.49391269E+00 

7 0.52955335E+00 

8 0.56020939E+00 

9 0.5903e647E+00 

0 0.62159282E+00 

1 0. 6523370 lE+00 


0.31250000E-01 
0.62500000E-01 
0.93750000E-01 
0. 12500000E+00 
0. 15&25000E+00 
0. 18750000E+00 
0,21875000E+00 
0.25000000E+00 
0.28125000E+00 
0,31250000E+00 
0.34375000E+00 
0.37500000E+00 
0.40625000E+00 
0.43750000E+00 
0.4&875000E+00 
0„50000000E+00 
0 . 53125000E + 00 
0.5fa250000E+00 
0.59375000E+00 
0 .62500000E + 00 
0.65625000E+00 


2 0.68312935E+00 

3 0,71397960E+00 

4 0.74490404E+00 

5 0.77591527E+00 

6 0.80703694E+00 

7 0. 83829 159E+00 
3 0.86972016E+00 
3 0. 90137768E+00 
} 0. 93336284E+00 
I 0.96588403E+00 
I 0 . 99999565E+00 


0.6875()0O0E + O0 
0.71875000E+00 
0 .75000000E+00 
0.78125000E+00 
0 .ei250000E+00 
0 .84375000E+00 
0 .87500000E+00 
0 . 9 0 6250 0 0 E + 0 0 
0 .93750000E+00 
0 .96S75000E+00 
0 . 99999899S+0O 





CHAPTER 3 


CLUTTER POWER SPECTRAL DENSITY AND SPATIAL CORRELATION 

ESTIMATION 


3.1 INTRODUCTION: 

In this chapter methods for the clutter power spectral 
density and spatial correlation estimation are presented. 
Adopting a scattering function model for a doubly spread target, 
an approximate model using tapped-delay lines is described. 

A method for finding the PSD of a clutter range cell using DpTs 
is given. The scattering function model assumes a sufficient 
spatial separation between the cells for their cross correlation 
to be negligible. As this need not be the case, two methods of 
estimating the cross correlation are also given. 

3.2 POWER SPECTRAL DENSITY OF THE CLUTTER FROM A RANGE CELL: 
ESTIMATION OF 

3.2.1 Scattering - function characterization of 

doubly spread targets : 

Clutter can be classified according to its spectral 
characteristics. The spectral width of distributed clutter 
is related to the nature of fluctuations of its component 
scatterers. A scattering function whose independent variables 
are both frequency (doppler shift) and range (or the two way 
time delay) can be derived starting from (1.1) [2], under the 
assumptions: 
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a) E[b(t, X )]= 0, where X is the two way time delay. (3.1) 

This assumption does not entail any loss of generality. 

b) The B processes from range intervals corresponding to delays A 

and are statistically uncorrelated 

E[b(t, X ) b*(u, Xj^)] = 0, Xt^X^ (3.2) 

c) The processes b (t, X ) and b(u,X 2 ) sire jointly WSS, 

We can define a two-variable covariance function K(Z,X ) 

E[b(t+Z, X) b*(t, Xjl)] = K(^,X ) ^(X-Xj^) (3.3) 

The covariance function of the clutter signal s(t) in (1.1) is 

OO OO 

// f(t^ E[b(tj^,Xj)b (tj.S)] 

•mOO 1.0 0 

f*(t2-X2)dXj^ dX2 (3.4) 

OO 

= Et / .f(t;^-X^)K (t;^-t2,Xj^)f (t2-*X^)dX3^ (3.5) 

wOO 

<• 

Thus, a second-order characterization of the process 
s(t) is obtained. 

The Fourier transform of K(Z,X) is the scattering 
function, 

S(f,X) =/ K(Z,X) dz (3.6) 

.. 1.00 

Some properties of the scattering function are: 


1. The average energy received is given by 



= / / / S(f,X) If(t-A) 1^ dfdXdt 


(3.8) 


whence 


E CO oo 

^ = / / S(f,A)dfdX (3.9) 

mmXX> mmOO 

since the energy of the transmitted signal is and 
oo _ 

/ lf(t)l2dt =1 (3.10) 

mmCO 

2* The normalised power density function in the f--)\ plane is 
given by 


9jq(f » ^) 


/ / S(f,X)dfdX 

wOO ..oo 


(3.11) 


Some models have been suggested in the literature for the 
clutter power spectral density 


OO 

p(f‘) = / gj^(f,X)dX 

■■■O O 

For example [7] 


P(f) 


1 

l+(f/f^)® 


where f is the upper cutoff frequency 
c 


(3.12) 


(3.13) 


a is a constant 


The value a = 3.0 was earlier suggested [8]. More recent 
work [ 7 ] indicates that a would vary with the type of clutter, 
and a of the order of a = 1.23 describes land clutter more 
closely than a = 3.0. 

3«2.2 The tapped-delay line model for a doubly-spread target: 

An approximate model for a doubly-spread target that is 
relatively easy to implement is a tapped-delay line model 
described here. 

The complex envelope of the composite clutter signal 
(1.1) is 

00 

s(t) = "fE. / f (t- A)b(t,X)d X , -<»<t<«' ( 3 . 14 ) 

«»co 

Assuming the transmitted signal f(t) is bandlimited, 

F{f^ = 0, If! > I (3.15) 

where F \ f j is the Fourier transform of f(t), we expand 
f(t-X) using the sampling theorem 


f(t-X) = E f(t - |-)( 
lc=— «) ”s 

where Wg = W 

Substituting (3.16) into (3.14) 


sin'rt.v^(A- k/Wg) 
itWg (X-k/Wg ■) 


) 


, we have 


s(t) = Z f(t - -|- ) b. (t) 

k=-~ Wg k 


(3.16) 


(3.17) 


where, 
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bj^(t)^ / b(t,X) 


simtWg(A-k/Wg) 

■ TtWgCX-k/Wg) 


(3.18) 


The target model given by (3.17) is shown in Fig. 3*1. 
The cross-corariance function of the tap gains, using (3.3) 
and (3.18), is given by 


^ * « simtvfg (X-k/Wg) simiWg(A- l/W^) 

E[b^(t) b^ (u)] =_/ K(t-u,X) (V-k7w ~ 1 • 

(3.19) 

Assuming K(t-u,X) to be constant with respect to X over l/W » 


E[bj^(t)bj_ (u)] = 


0, k 4: 

ijj|^ K(t-u, k/W ) , k = 4t 


(3.20) 


The tap gain processes can be characterixed in terms of 
their spectra, by taking the Fourier transform of (3.20), 

S|j(f) = V V ) 

s s 

Sj^(f) is just the cross-section of the scattering function 
S(f,X) at various values of X . 

We, thus, have an approximate model for the target, 

R 

?p^(t) = I f(t - ^ )b. (t) (3.22) 

k=o s 

3.2.3 Estimation of relative power spectral density from the 
data: 

The power spectral estimates of the clutter data measured 






from the Kth range cell are obtained by finding the square of 
the modulus of the DFT of the bjj,(n)^ sequence, i*e* 


n = 1,«.«,N m = (3*23) 

Gj^(m) = lBj,(m)l^, m = (3.24) 


where Gj^(m) denotes the relative power spectrum. 

It is convenient to write (3.24) in the form of a 
probability density 

- _ 

I lBj^(m)l^ 

m = o ^ 



(3.25) 


The cumulative density function can be similarly found. 
The various steps in the evaluation of the normalised PSD 
function are shown as a flow chart in Fig. 3,2. 

A radix 2, in-place, decimation-in-time FFT algorithm has 
been used for the computation of the DFT and is shown in the form 
of a flow-chart in Fig. 3.3. The algorithm is explained in 
terms of a decimation-in-time signal flow graph in Fig. 3.4. 

The resolution in both frequency and range, with respect 
to the measurement system in Table 1.1, obtained are: 

a) range resolution s ^ T^ 

where C is the velocity of light 

Tq is the range interval in two-way time delay units 




Fig. 3.2: EVALUATION OF IHE POWER SPECTRAL DENSITY OF THE 
CLUTTER FROM ONE RANGE CELL. 








Fig. 3.3: FFT ALGORITHM 
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FI6.3.4 SIGNAL FLOW GRAPH (DECIMATION - IN -TIME) FOR THE FFT(N=16) 
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For = 0.1 jis, = 15 m 

b) frequency resolution ® ^ 

where T is the interval between samples, T = 3,27670 sec 
N is the no. of samples in the sequence, N = 2^,r 
being an integer. 

The highest recoverable frequency fpj = ^ 

= 152.59 Hz. 

For N = 1024, f^ = 0.298 Hz 

An illustrative example of PSD estimation using a 
program incorporating the scheme of Fig. 3,2 is shown at 
the end of the chapter, 

3.3 SPATIAL CORRELATION ESTIMATION; 

3,3.1 Cross Correlation Function; 


The cross correlation between the clutter from any two 
range cells is given by the function 

\,L^^1»^2^ = E [bj,(tj^). bj^ (^2^^ (3.26) 

where K and L are the range cell indices. 

Assuming the processes to be jointly wide-sense stationary, 
we can write (3.26) as 

= E[bj^ (t + u) bj^ (t)] 

In general Rj^ would depend upon: 


(3.27) 



1. the relative rate of fluctuation, density and nature of the 
point scatterers in the two range cells, 

2. the spatial separation between the cells, 

3. the temporal separation u. 

A study of the correlation properties of the various sources 
of clutter is useful in the estimation of clutter from a parti- 
cular range cell by using the clutter data obtained from adjoi- 
ning cell(s). These estimates could be used in the detection of 
weak radar signals, for the purpose of clutter cancellation in 
the cell of interest, 

3.3.2 Evaluation by the Direct Method: 

The normalised cross covariance (called correlation 
coefficient henceforth, after Papoulis [6]) between the 
continuous complex random processes b^, and bj^ is the quantity 
to be evaluated and is given by 


‘ {E[|b^(t+u)-i:^l2].E[|bL(t)-5Ll^]l^^^ 

where, 

Pj^(t+u) = E[bj^(t+u)] 

= E[bj^(t) ] 


The data consists of the two sequences 

f** 

Kth range cell : bj^(n), n = 1,,..,N 
Lth range cell : bj^(n), n = 1,...,N 


(3.28) 


(3.29) 
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The value of 5 is estimated as 


J CM) _ n=o 
K L ^ 

|[ |b^(n+Nj^)-tl^|2][ lbL(n+N;^+M)-pLl^]] 


( 3 , 30 ) 


where. 


jijj, = Y 2 bjjrCn+N^) 


= Y ^ b^Cn+Nj^+M) 
n=o 

Nj^,M and L are variable parameters having integer values: 

Nj^ : index of the first sample to be taken in the first 
sequence, for the computation of S » S. ^ 

M : no* of samples corresponding to the time difference 
u = ° i ^ 1 N-Nj^. 

I : no. of samples corresponding to the time frame over 
which I is to be evaluated, 2_<L^N-Nj^-M+1 

Varying N^, in effect, enables a shift along the time- 
axis, of the interval over which J is estimated, while varying 
I varies the length of this interval. Varying M, in effect, 
changes the time difference (u = t^^-t^) in increments or 
decrements equal to integer multiples of the sample— to— sample 
interval. 
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Fig. 3.5: EVALUATION OF CROSS COVARIANCE (NORMALISED) BETWEEN 
THE COMPLEX REFLECTIVITY PROCESS OF TWO RANGE CELLS: 
DIRECT METHOD 
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The various steps in the evaluation of 
are shown in the flow chart of Fig. 3.5. 

3.3.3 Evaluation by the FpT Method: 

The direct (conventional) method described in the 
previous section is suitable when the no. of data points used in 
the computation (I) is small, i.e. when I < 64, as a thumb-rule [ 10 ]< 
Again, only one value of ^ computed at a time» corres- 

ponding to a preselected value of M and I. In case a series 
of values 

is required, the method is cumbersome as Nj^, M and I have to 
be fed in each time. 

When it is required to evaluate the correlation function 
over all possible values of M and for large values of I (I > 64. 
as a thumb-rule), a faster method employing Fast Fourier 
Transforms can be used. 


For convenience we change the symbols for the two 
sequences as follows, and we specify N = 2“^ , where J can take 
non-negative integer values. 


bj^(n) = x(n) 
bL(n) = y(n) 


n = 1,...,N; X and y denoting 
complex numbers 


(3.31) 


j 

Since sample indices are numbered n=l,...,N, rather than 
from n=0, the Discrete Fourier Transform (DFT) can be written as 
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F[x(n)] = E x(n) l<k<N (3.32) 

n=l 

l<n<N 

4 X(k), 1 < k < N 

The correlation function corresponding to the changed 
notation (3.31) is 

1 N ♦ 

z(i) = w E x(n+ 4 ) y (n), 0 < I < N-1 (3.33) 

n=l “ " 

The DFT pair corresponding to discrete-time correlation 
is 

i E x(n+4) y*(n)<^=>i X(k) Y*(k) (3.34) 

n=l 


0<1<N-1, l<k<N 

1 < n < N 


In order to use (3.34) the two data sequences are assumed 
cyclic with period N, i.e. 


x(n+N) = x(n) 
y(n+N) = y(n) 



1 f 


,N 


(3.35) 


The normalised value of the cross-correlation 


'( 4 ) = 


m. 


1 

N 


N o N 

[ E !x(n)l^][ E ly(n) 
n--=l n=l 


1/2 


, 0 < i < N-1 


(3o36) 


is evaluated here (as compared to the normalised crosscovariance 
in the first method) through the steps shown in Fig. 3.6. The 
specific FFT algorithm used is the same as for PSD estimation. 




READ 

X ( n ) — b|^ 

y(n) = Bj^(n), n=li...,N 


OTPUTE ; : 

X(k) = F[x(n)],lc=l,...,N 

Y(k) = F[y(n)],k=U...,N 

1 

1 — : 

COMPUTE 

Dj^ = Z lX(k)l2 

D2 = T |Y(k)l^ 

^ k=l 

' — . , , ' ,4^ '... , — . 

D = V'Djl.D2 • 

^ ~ 

COMPUTE 

W(k)=X*(k) .Y(k) , k=l,...,N 
COMPUTE P(l)= F[Ki(k)] *, 1=0,..., N-1 
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NOTE: 1. The DFT and the inverse DFT are both computed by the same 
’fast’ algorithm by meafuSs* of the relation 


F"^[x(k)]= i[ 


k=l 


F[X*(k)] * , n=l,...,N 


»N 


2. x(n) = bj^(n) 
y(n) = bj^(n) 


n=l , . . . ,N 


LEGEND: F[.] denotes the DFT of the sequence within brackets. 

Fig. 3.6: Evaluation of cross-correlation (Normalized) between the 
complex reflectivity processess of two range cells: 

FFT method. 
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3*4 SOME ILLUSTRATIVE EXAMPLES USING GENERATED DATA 
3.4.1 Estimation of correlation from the datai 
The following cases are considered: 

1. Autocorrelation of the data from one range cell (by th'e 
EFT method). 

2. Cross correlation between the data from two range cells 

a) when positively perfectly correlated 

b) when negatively perfectly correlated 
(by both methods). 

3. Cross correlation between the data from two range cells; when 
the value of correlation is low (by the FFT method) . 

Example 1: 

The autocorrelation properties of the data obtained 
from a particular range can be investigated. We consider the 
data of example 2, Chapter II (Table 2.5), where the 
magnitude values follow a lognormal distribution. The auto- 
correlation of this complex-valued sequence, evaluated by 
the FFT method, is shown in Table 3.1. As expected, normalised 
correlation values are symmetrical about the value of time 
shift index 

M = ™ = 16 
2 

Since the methods assumes the sequence to be periodic. 

‘ ' ' ' 

* ! 

Acc. No. 
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TABLE 3.1 

CORE E L A I I 0 N ( F F T MET H 0 B 


r is the c o r r e 1 a t i o n v a i u e ( normalised ) . 


B e < r ) 


0.99gg9Aq4E+00 
0.10448679E+00 
0 . 17978633E-01 
-0.10911529E-01 
0.3032851SE-01 
0.901S&871E-02 
-0 . 19834511E-01 
0.2&362482E-03 
0.15620577E-01 
-0.16&32585E-01 
0.11454485E-02 
-0.206067S5E-02 
0.S3971350E-02 
0,80722850E-03 
0.93083074E-02 
0.19427773E-01 
0.22854383E-01 
0.19427321E-01 
0.9308305GE-02 
. 0.80722431E-03 
0.83971424E-02 
-0.20&0&755E-02 
0.11454439E-02 
-0.16&32600E-01 
0. 15&20572E-01 
0.26363973E-03 
-0.19834507E-01 
0.90166861E-02 
0.30328525E-0i 
-0. 10911531E-01 
0. 17973&38E-01 
0. 1044868SE+00 


I m ( r ) 


O.OOOOOOOOE+00 
0.22927677E-01 
0.43906569E-01 
-0.23324009E-01 
0-52gi2477E-02 
0.20706305E-01 
-0.595G4S73E~02 
-0. 19726947E-01 
0.184565918-01 
0. 13054a60E-01 
0.1i694972E~01 
-0.68387045E-02 
0.45178696E-02 
-0. 122526fa2E-01 
-0. 12592740E-01 
-0.54777153E-02 
O.OOOOOOOOE+OO 
0.54777395E-02 
0. i2592763E-01 
0. 12252b44E-0i 
-0.45178S2.6E-02 
0.683S7082E-02 
-0. 11694972E-01 
-0. 13054S51E-01 
-0.18456591E-01 
0.1972fa932E-01 
0 .59564891E-02 
-0 .20706305E-01 
-0.52912547E-02 
0,23324024E-01 
-0.43g0659bE-01 
-0 .22927696E-01 


Abs ( r ) 


0. 999999g4E+00 
0. 106972758+00 
0.47444895E-01 
0.257501&2E-01 
0.30786626E-01 
0.22584327E-01 
0. 207096028-01 
0. 1972S709E-01 
0 .24,179500E-0i 
0.21144083E-01 
0. 11750934E-01 
0.7142426gE“02 
0.95353564E-02 
0. 12279224E-01 
0 . 15659556E-01 
0 .201S5236E-01 
0.228543S8E-01 
0 ,20185290E-01 
0.15659573E-01 
0 .12279205E-01 
0.95353592E-02 
0 .7i42430&E-02 
0.11750932E-01 
0.21144090E-01 
0.24179496E-01 
0.ig728694E-01 
0.20709598E-01 
0.22584327E-01 
0.30786633E-01 
0.25750175E-01 
0.47444925E-01 
0,10697283E+00 


Phase(r) (deq) 


0 .OOOOOOOOE + OO 
0.12376337E+02 
0.67732132E+G2 
0.64928680E+02 
0 .98964691E+01 
0.&6468987E+02 
0.16715500E+02 
-0,8g234360E+02 
0.49757378E+02 
-0.38i28246E+02 
0.84406082E+02 
0.73231102E+02 
0. 2828136 8 E+ 02 
-0.86230690E+02 
-0.53528912E+03 
-0. 15745947E+02 
0. OOOOOOOOE+OO 
0. 15745977E+02 
0.53528969E+02 
0.S6230705E+02 
-0.28281311E+02 
-0.73231117E+02 
-0. 84406097E+02 
0,3B128204E+02 
-0.49757389E+02 
0.S9234314E+02 
-0.1&71550SE+02 
-0.66463987E+02 
~0.gsg64787E+01 
-0.64928688E+02 
-0. 67732 140E+02 
-0.12376337E+02 


Re(.), Im(.), Abs(.) denote the real part, imaginary part 
and absolute value respectively. 
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Example 2: 

We consider a hypothetical case where the clutter 
variations as measured from two range cells K and L are 
similar in all respects, except, due to the relative locations 
of K and L with respect to the radar, for a scaling factor 
S = 3.0. The data from cell K is taken as the data of 
example 1, Chapter II (Table 2.2), where the magnitude values 
follow a ^ayleigh distribution. The number of samples (N)=32. 

The data from cell L is shown in Table 3.2. ftesults of 
correlation performed by the direct method, for a few values 
of the parameters N^, M and L, are shown in Table 3.3, 
where, 

= starting sample number in the first sequence; 

1< Ni < N 

M = time delay; 0 <_ M N - 1 

L = number of consecutive samples over which correlation 
is evaluated; 1<.L^N, 

subject to the condition + M + L-1 < N, which is checked 

by the program, f^esults of correlation performed by the FFT 

method are shown in Table 3.4. For M 0, there is a difference 

as expected in the values given by the two methods, since 

normalised cross covariance and normalised cross correlation are 

respectively arrived at by the direct and FFT program segments, and 
since dc content is present in the input sequence. (In the FFT 

program segment, Nj^=l and L=N are fixed values, and correlation 

is evaluated for all M, 0 _< M £ N-1). 
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TABLE 3.2 
CELL DATA 


IN PHASE PART 

1 0.83735199E^-01 

2 0. lb593754E+02 

3 0.&3011022E+01 

4 ■ 0. 14702725E+01 

5 -0,27790837E+01 

6 0.&175732SE+01 

7 0.24207370E+01 

8 ~0.5292720aE4-01 

9 -0,18921965E+00 

10 0 .38407242E+01 

11 -0.39480035E+01 

12 -0.11367635E+01 

13 0.43451491E+00 

14 0.11932796E+01 

15 -0.14746627E+01 

16 0.40&98533E+01 

17 0.35722437E+01 

18 0.36613345E+01 

19 0.24810073E+01 

20 -0.2001232&E+01 

21 0.2653S374E+01 

22 0.38816822E+00 

23 -0.50204134E+00 

24 0.10025810E+01 

25 -0.17138822E+01 

26 -0.42642385E+00 

27 0.17295971E+01 

28 0.17970763E+01 

29 -0.14098&59E+01 

30 -0. 15038581E+01 

31 0.74683S09E+00 

32 ' 0.94243771E+00 


QUADRATURE PART 
0 . 25563064E+00 
0 . 20357499E + 0 i 
0.28688633E+01 
0 . 73473282E+01 
-0.54650474E+01 
0. 19922533E+01 
0.49916439E+01 
-0.24259326E+01 
-0.50647368E+01 
0 . 3649S349E+01 
0.2457038&E+01 
0.471S2178E+01 
-0.42487903E+01 
0.42939358E+01 
-0 .3628S435E+01 
-0.42015955E+00 
0 .218308S1E+00 
0.7942696SE+00 
0.2097824eE+01 
0.27651749EH-01 
0 . 12193145E + 01 
0 .30605390E + 01 
-0 .25367994E+01 
0.2565&142E+01 
-0.143580i6E+01 
-0.23756S30E+01 
0.67003828E+00 
-O . 98752940E + 00 
-0.75428672E-01 
0.6651486SE+00 
0.29i79499E+00 
0„64708340E+00 
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TABLE 3.3 

Cross-covariance (Normalised) (r) by Direct Method (Example 2) 


(a) M = 0, = 1 


L 

Re(r) 

Im(r) 

Abs(r) 

Ph (r) (deg) 

8 

O.IOOCOOOOE+Ol 

-0.38345824E-08 

O.lOOOOOOOE+01 

-0.21970537E-06 

.6 

O.lOOOOOOOE+01 

0.60278986E-08 

O.lOOOOOOOE+01 

0.34537314E-06 

»4 

0.10000002E+01 

-0.15261712E-10 

0.10000002^+01 

-0. 87443142 E-09 

32 

O.lOOOOOOOE+01 

0.65132508E-08 

O.lOOOOOOOE+01 

0.37318176E-06 




(b) = 

1, L = 16 


M 

^e(r) 

Im(r) 

Abs (r) 

Ph(r) (deg) 

4 

0.27008435E+00 

-0.17715546E+00 

0.32300094E+00 

-0.33261925E+02 

8 

0.21851349E+00 

-0.28608727E+00 

0.35999176E+00 

-0. 52627403 E+02 

12 

0.56510989E-01 

-0.10714775E-01 

0.57517804E-01 

-0.107361 31 E+02 

16 

0.22118418E+00 

-0.58230016E-01 

0.22 872075 E+00 

-0. 14749292^+02 



(c) M = 0 

, L = 16 


M| 

Re (r) 

Im(r) 

Abs (r) 

Ph(r) (deg) 

2 

0.99999994E+00 

0.23430551E-08 

O.99999994E+OO 

0.12851760E-06 

4 

0. 10000002 E+01 

-0.18534123E-08 

0.10000002 E+01 

-0.1061 92 67 E-06 

8 

O.99999994R+OO 

0,42214330E-08 

O.99999994E+OO 

0.24187031E-06 

L 6 

0.10000001 E+01 

-0.63362982E-08 

O.lOOOOOOlE+01 

-0.36304306E-06 


'lote: Re(.), Im(.), Abs ( . ) and Ph(.) denote the real part, imaginary 
part, absolute value and phase respectively. 
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TABLE 3„4 

C Q E: E: E L A T I ON! ( F E T 


r is t h e 
Ee ( r ) 


O.lOOOOOOOE+01 
0 .2121S443E+00 
0.30918017E-01 
0 .96S74092E-01 
0.15747385E+00 
-0.64&74772E-01 
-0.55691&82E-01 
0 ,314S94.24E-01 
0.i0393536E";-00 
“O > 17b72403E + 00 
0,64104125E-01 
0.760618SOE-02 
0-11210383E+00 
-0.42503420E-01 
0.1782S340E+00 
0„283&07i8E-K)0 
0.20140755E^-00 
0.28360704E+00 
0. 17828S33E+00 
-0.42503323E-01 
0.11210384E+00 
0.760&1189E-02 
0.64104088E-01 
-0. 1757240:3E + 00 
0. 10393530E+00 
0.31489499E-01 
-0.55&91b71E-01 
-0.&4674839E-01 
0. 15747385E+00 
0.96t.742igE-01 
0.30917995E-G1 
0.21218443E+00 


correlation value 
I m C r ) 


-0.63r333991E-0S 
0.316744f30E + 00 
0.38S81302E-01 
-0. 12555692E+00 
0„912fa4315E~01 
0. 16&92807E+00 
-0.94160937E-01 
- 0 . r3t309625E-^-00 
0.21296953E+00 
0.5237S147E-01 
0.342635928-01 
0.20459294E-02 
~0.28170757E-01 
-0. 10358461E+00 
-0. 19614510E+00 
-0.33399165E-01 
0.61459562E-0S 
0.33399343E~01 
0. 19614497E+00 
0. 10358448E+00 
0.28170794E-01 
-0.20459145E-02 
-0.34263633E-01 
-0.52376114E-01 
-0.21296950E+00 
0.i3609622E+00 
0 .94160996E-01 
-0 . 16692810E + 00 
-0.91264337E-01 
0 . 12555696E+00 
-0 .38381 190E-01 
-0 .31674477E+00 


C ri 0 r m a 1 i s e d ) 
Abs(r ) 


0. lOOOOOOOE+01 
0,33124713E+00 
0.49675744E-01 
0. 15846267E+00 
0. 18200876E+00 
0. 17901900E+00 
0. 10939765E+00 
0.13969171E+00 
0 „ 23 6 97 8 02 E+ 00 
0. 18336357E+00 
0 .72686538E-01 
0 .78765405E-02 
0 . 11558920E+00 
0. 11 196567E+00 
0-26506537E+00 
0.28556705E+00 
0 .20140755E+00 
0.28556693E+00 
0.26506522E+00 
0.11196552E+00 
0.11558921E+00 
0,78764716E-02 
0.726S6531E-01 
0.18336356E+00 
0.23697796E+00 
0.13969171E+00 
0.10939769E+00 
0.17901906E+00 
0. 18200877E+00 
0. 15846279E4-00 
0.49&75643E-01 
0.38124728E4-00 


Phase ( r ) ( deg ) 


-0 .36459591E-06 
0.56182217E+02 
0.,51508636E+02 
-0.52405083E+02 
0.30094528E+02 
-0.&S821587E+02 
0.59397720E+02 
“0.76972366E+02 
0 .63986267H+02 
-0.1 6597 lfaSE+02 
0.28124462E+02 
0. 15055211E+02 
-0. 14105874E+02 
0.67690392E+02 
-0.47730358E+02 
~0.67165351E+01 
0. 17433S20E~05 
0.67165737E+01 
0.47730354EH-02 
-0.67690414E+02 
0. 14105891E4-02 
-0. 15055234E+02 
-0.28124502E+02 
0. 16597157E+02 
-•0.&3986275E + 02 
0.76972336E+02 
-0.59397739E+02 
0„ 68821 571 E+02 
-0.30094536E+02 
0.52405C56E-i-02 
-0.51508575E4-02 

— U M 5 w 1 8 2 3 11 1 V 2 
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Example 3: 

Here we consider two range cells 0 and P whose clutter 
returns are 180® out of phase but otherwise similar. The data 
from Cell 0 is shown in Table 3.5. The magnitude values of 
this cell data are Weibull distributed with parameters = 1,0 
and a = 3,0. The data from cell P, shown in Table 3.6, is a 
scaled by a factor of S = --0,2 version of the same data. 

Results of correlation evaluated by the two methods are shown 
at 

(a) table 3.7 for the direct method (for various M and L 

values ) 

(b) table 3.8 for the FFT method (for all values of M, with 

= 1 and L = N). 

Example 4: 

While examples 2 and 3 were hypothetical cases meant 
primarily to check the program functioning, here we consider a 
case more likely to represent a real situation. We consider 
two range cells Q and R, The data consists of N = 32 complex 
samples. The magnitude values follow a Weibull distribution, 
but with different parameter values for the two cells: 

Cell Q : p = 1,0 o = 3.0 (Table 3,5, example 3) 

Cell R : p = 4.0 o = 7.0 (Table 3.9) 

The FFT method is used to obtain estimates of the normalised cross 
correlation, and the values obtained are given in Table 3.10, 



TABLE 3.5 
CELL DATA 

IN-PHASE PART QUADRATURE PART 
1 0. 10332370E-f02 


2 0 . 41038S35E+02 

3 0.64630241E+01 

4 0. lfa321096E+01 

5 -0.25242&03E+01 

3 0.53370475E+01 

7 0 . 19895430E+01 

8 ~0.45&52270E+01 

9 -0 . 14207643E+00 

10 0,3014738SE+01 

11 - 0.271 98 l&2E+0i 

12 -0 .817329G5E+00 

13 0 .27493137E i-00 

14 0.787S5747E+00 

15 -0.85574S41E-f00 

13 0.246S9244EH-01 

17 0. 18g40344E+0i 

18 0.20321739E+01 

19 0. 1 1942064E-f01 

20 -0.1011993&E+01 

21 0. il482440E+01 

22 0. 17741053E+-00 

23 -0. 19233757E+00 

24 0. 40913469E+00 

25 -0. 587fa9508E+00 

26 -0. 15247783E+00 

27 0.4752793gE+00 

28 0 . 54592150E+00 

29 -0.29489833E+00 

30 — 0 . 3 b 6 3 5 9 4 7 E + 0 0 

31 0.S8715233E-0i 

32 0. 15961395E+00 


0.31726301E+00 
0.50420752E+01 
0.29425857E+01 
0 .815&0698E+01 
-0 .49639392E+01 
0. 19152551E+01 
0.41025066E+01 
-0 .20924840E+01 
-0 .38028808E+01 
0 .28&49023E+01 
0.16g26768E+01 
0.33923824E+01 
-0.2bS83445E+01. 

0.283505i8E+01 
-0.2105S221E+01 
-0.25467795E+00 
0.11574912E+00 
0. 44084862E+00 
0. 10097655E-i-Oi 
0. 13983078E+01 
0.52756453E+00 
0 . 13938055E+01 
-0.97187585E+00 
0. 10469795E+01 
-0.47558549E+00 
-O.S4947026E+00 
0.1S4121i5E+00 
-0 .29999480E+00 
-0,15777240E-01 
0,16203891E+00 
0.34661b73E-01 
0.1095918gE+00 
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TABLE S.ri 
CELL DA 


IN-PHASE PART 

1 -0.20784740E+01 

2 -0 .8219767SE+01 

3 -0.1292&048E+01 

4 -0 . 32642195E+00 

5 0.5048520&E+00 

6 -0.11874095E+01 

7 -0.397908&OE+00 

8 0.91304541E+00 

9 0.284152S7E-01 

10 -0.S0294777E+00 

11 0.5439S325E+00 

12 0. 1634658 lE+00 

13 -0.54986276E-01 

14 -0. 15757i49E+00 

15 0. 17114969E+00 
IS - 0 . 4 9 3 3 8 4 9 0 E + 0 0 

17 - 0 , 3 7 8 8 0 6 8 9 E + 0 0 

18 -0.40643477E+00 

19 -0.23884128E+00 

20 0.20239873E+00 

21 -0.22964881E+00 

22 -0. 354821 05E-01 

23 0. 384675 15E-01 

24 -0.8182S940E-01 

25 0.11353902E+00 

26 0.30495567E-01 

27 -0.95055878E-01 

28 -0.10918430E+00 

29 0 .58979668E-01 

30 0.73271893E-01 

31 -0 . 17743047E-01 

32 -0.31922791E-01 


QUADRATURE PART 
-0.&3452601E-01 
-0 . 100S4151E+01 
-0.58851713E+00 
-0. i6312140E+01 
0.9927S7S4E+00 
-0.3S305101E+00 
-0 .82050133E+00 
0 .41849682E+00 
0 .76057619E+00 
-0.57298046E+00 
-0.33S53537E+00 
-0. 6784765 lE+00 
0.53766888E+00 
-0.56701034E+00 
0 . 421 16442E+00 
0.50935593E-01 
-0.23i49824E-01 
-C.8S169724E-01 
-0 , 20195310E+00 
-0.27966157E+00 

— 0 . 1 0 U 5 1 2 9 1 E + 0 0 
-0. 279761 llE+00 

0 . 19437517E + 00 
-0 .20939592E+00 
0,951 171 OOE-01 
0.169S9405E+00 

— 0 .36824230E— 01 
0,5999S959E-01 
0,315544S0E-02 

-0.32407783E-01 

-0,69323345E-02 

-0,21918377E~01 
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TABLE 3.7 

ross-covariance (Normalised) (r) by Direct Method (Example 3) 


L Ee(r) Im(r) Abs(r) Ph(r)’(deg) 


D 32 
L 31 
L 30 
I 29 
> 28 
} 27 
I 26 
^ 25 


-0.1 0000001 E+01 
-0.26604116E+00 
-0.24358593E+00 
-0.57376287E-03 
0.32842219E+00 
-0.38023984E+00 
-0.26693031E+00 
-0.14688352E+00 


-0.10531747E-07 
•0.14 80411 9 E+00 
0.30391911E+00 
0.30231792 E+00 
-0.32444155E+00 
-0.31417581 E+00 
-0. 35366309^+00 
0.33424926E+00 


0. 10000001 E+01 
0.30445704E+00 
0.3 8 948807 E+00 
0.30231848E+00 
0.46165296E+00 
0.49324310E+00 
0.44309074E+00 
0.36509910E+00 


0.60342455E-06 
-0.29094175E+02 
-0.51288425 £+02 
-0.89891251E+02 
-0.44650658E+02 
0.39565456E|02 
0.52955997E+02 
-0. 66277245 E+02 


Note: Re(.), Im(.), Abs(.) and Ph(.) denote the real part, 
imaginary part, absolute value and phase respectively 
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)< S e ( r ) 

I m ( r ) 

A b s ( r ) 

PheseCr ) (de*^) 

-0 . lOOOOOOlE+01 

-0..37221760E-03 

0. lOOOOOOlE+Ol 

0 .55703958E-06 

-0 .307513&3E+00 

-0.i3710643E+0Q 

0»33669403E-^*00 

0 . 24029934E+02 

-0 .4S343Sif3E-01 

-0 . 12228b04E + 00 

0. 13167989E-i-00 

0 .&8227112H-^02 

0 .20435425E-02 

0„859i3867E-01 

0.85933171E-01 

0 .8863742 lS-f-02 

-0 . 12307041E+00 

-Q .34315757E-01 

0. 12776501E+00 

0 « 155800793-^02 

-0.22129023E-01 

— A Q *7 I i *4 fZ t; A 7 

w • -f / i. .L ’7; w O Zt s/ i. 

0.99603996E-01 

0.77163475E+02 

0.6808.1960E-01 

0.3&217906E-01 

0.77116087E-01 

0.2301 1864E+02 

0 .58417339E-02 

0.77528045E-Q1 

0.77747822E-01 

0.S5690910E+O2 

-0 .&915r975E-01 

-0.79673833E~01 

0, 10550220E+00 

0 .49045315E+02 

0. 657386 18E-01 

-0.533578S1E--01 

0.84&&7757E-0i 

-0 .390&5105E't-02 

-0.6ir95716E-02 

-0.48783273E-01 

0.4916560bE-01 

0 - 3284 93 3 OH + 02 

0.63128937E-02 

0. 19574270S-01 

0.205S7078E-01 

0 « 72124377E-5-02 

-0 .41378930E-01 

-0=37245434E-02 

0.41 5462 14E-01 

0.51433682K+01 

0 . 18722257E--02 

0 „ 46327143E-0i 

0.4&364.959E-01 

0 .87 6 8 5 7 5 3 E + 0 2 

-0.52954<)89E-01 

0.66271946E-01 

0.S4829867E-01 

-0.51373665E-f02 

-0.39242724E-01 

0.24040259E-01 

0.10211293E+00 

-0. 13616860E+02 

-0.10419881E+00 

-0. 13445560E--07 

0.10419381E-fOO 

0.73933070E--05 

-0.99242792E-01 

-0.24040390E-01 

0 .10211304E+00 

0. 136i6923H+02 

-0.5295404aE~01 

-0. G6271976E-01 

0.84829367E-01 

0.51373695E-^02 

0„lS722039E-62 

-0.46327055E-01 

0 .46364866E-01 

-0.8763577SE+02 

-0.41378967E-01 

0,37245322E~02 

0.41546252E-01 

-0.51433437E+01 

0.S3128i353E-02 

-Ou 1957426SE-01 

0 .20567065E-01 

-0.72125053E-5-02 

-0.61195884E-02 

0. 48783258E-01 

0 ,49165592E-01 

-0. 32349907E+02 

0.657386S3E-01 

0,53357873E-01 

0 .34667794E-01 

0.39065079E+02 

-0.69151931E-ai 

0,.79fa7B833E-01 

0.105S0217E+00 

-0. 49045334E+02 

0„58416724H-02 

-0,. 77527978E-01 

0 .77747747E-01 

-0 . 356909488+02 

0.6S081953E-01 

-0.36217898E-01 

0.77116080E-01 

-0. 230118 ■-02 

-0.221290I0E-01 

0.97114697E-01 

0.9960400:3E-01 

~ 0.77163-^’ 02 

-0. 12307042E+00 

0.34315805E-01 

0.12776501E+00 

-0.15580. ;j--02 

0.20435220E-02 

-0„85913897E-01 

0-8593S200E-01 

-0.3S637436E+02 

-0. 48843820E--01 

0. 12228603E+00 

0.13167992E^00 

-0. 682271 19E+02 

-0.30751389E+00 

0. 13710658E+00 

0.33669424H+00 

-0.24029933Ef02 
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r is the correlation value < normalised ). 


Re ( r ) 


-0 « 32323885E-i-00 
-0.11239727E+00 
-0.256974&iE+00 
-0.3280390SE+00 
0. 10373923E'^-00 
0.2054&307E+00 
-0. 1<3173336E + 00 
-0.25239691E+00 
0.59156001E-01 
0.2315g978E+00 
-0.40732g35E-01 
-0.32184020E-01 
-0. 10011426E+00 
-0.30579131E-01 
-0.85467711E-01 
-0.22311495E + 0,0 
-0. 17G08058E+00 
-0.28238043E+00 
-0. 13290328E+Q0 
0.51903&05E-01 
-0.78484848E-01 
0. 19362621E-01 
-0.b4399838E-01 
0. 15022567E + 00 
-0 . 12759882E+00 
0. 145bl258E-01 
0.11202595E+00 
-0 .9482&9f37E“02 
-0.131029a3E+00 
-0 .27346574E-02 
-0 . 1054S800E-01 
-0 .&9952682E-01 


I m < r ) 


-0. 11175S71E-06 
-0.770545SOE-01 
0. 13394741E+0Q 
0,11193673E+00 
0. 40536571E-01 
-O.49415704E-0i 
0. 13485540E-0i 
0.92034005E-02 
-0.29147747E+00 
-0,44S77812E-01 
0.i2412541E+00 
-0.151936&2E+00 
0.84999911E-01 
0.i048S811E-i-00 
0 .25320908E+00 
0.55078831E-01 
0.35307392E-01 
-0.27339198E-02 
-0 .86243965E-01 
-0 .il271880E + 00 
0 .50584171E--01 
-0 .9862&427E-01 
0.95342i07E-01 
0,&7049429E-01 
0 .12549472E + 00 
-0.17772b&4E+00 
-0,591&S439E-01 
0.13976103E+00 
0.55986207E-01 
-0.10144290E+00 
0.50333917E-01 
0.132S18&1E+00 


A b s ( r ) 


0 .323238808"^ U 0 
0 . 13!327380E + 00 
0.29213411E+00 
0. 34661 141E+00 
0.11604927E-4-00 
0.21132i99E+00 
0 . 16229450E+00 
0.25256464E+00 
0.29741932E+00 
0-23590778E+00 
0.13063802E+00 
0.15530790E+00 
0.131331068+00 
0,10g25475E+00 
0.26724443E+00 
0.234S7019E+00 
0. 17958553E+00 
0 . 2S239366E + 00 
0,20222132E+00 
0. 12409476E+00 
0 . 93373604E--01 
0. 10050gilE + 00 
0. 11505415E+00 
0.. 16450g51E + 00 
0. 17397034E+00 
0. 17832215E+00 
0. 126S9049E+00 
0. 14003234E+00 
0. 1424S955E+00 
0. 10147g75E + 00 
0.51427018E-01 
0. 149g3bg0E + 00 


Phase(r) (deq) 


0 It 1. w 0 ^ 0 

0.34432800E+02 
-0 . 28400305E+02 
-0. 18841 173E+02 
0. 20444301 E+02 
-0. 13523302E+02 
-0.47&63713E+01 
-0.20883081E+01 
-0. 78527512E+02 
“ U > 1 0 g 6 6 4 7 4 E + 0 2 
-0.7i832245E+02 
0.7S040092E+02 
-0.40332222E+02 
-0.73746460E+02 
-0.71348503E+02 
-0.1 3 574 373 E+02 
-0. 1133S488E+02 
0 .55470258E. + 00 
0.25244""' 02 

-0.65275 "02 

-0.32 302 ..;2£+02 
-0 .7S3g2'^g2E + 02 
-0.559625242+02 
0.24052387H+02 
-0 .4452367e''' + 02 
-0.853161^" -02 

-0.2784C .+02 

-0.86111 .E+02 

-0.23135 j9E+02 

0.88455313E+02 
-0.73165649E+02 
-0.6218g579E+02 
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As stated earlier the sequences are assumed to be periodic 
(with periodicity N) and correlation is performed over N samples* 
The results show a low value of correlation between the two 
sequences, 

3,4.2 Estimation of PSD from the data: 

In the following example, the complex reflectivity data 
consists of a relatively slowly varying sequence, whose PSD is 
similar to a low-pass-filter transfer characteristic. 

Example : 

The data sequence considered (Table 3.11) consists of 
N = 32 samples of the complex reflectivity of a given range cell, 
in which successive sample magnitudes decay exponentially 
over the interval 0 <_ t) £ N/2. The data is symmetric around 
n = ^ = 16. Normalised values of the PSD obtained with the 
program are shown at Table 3.12. The PSD values are plotted in 
Fig, 3.7 for the frequency index 0 <_ k £ 16, where ^ = 16 
represents the maximum recoverable frequency. The cell data 
are also plotted in Fig. 3.7. 

Normalised PSD values were also obtained, for purpose of 
comparison, through direct computation of the DFT (i.e. by 'brute 
force' method) and these are shown in Table 3,13. As expected, 
it is seen that the two sets of values of normalised power 
spectral density in ^ables 3.12 and 3.13 are equivalent. 



TABLE 


3.11 


CELL DATA '■ 

X ( rr ) = R e ( K < n ) ) + ,j A I m < x ( n ) ) 
Abs(x(ri)) is its magnitude 
The complex sequence [K(n)j represan 
varying clutter signal ; 
n Re(x(ri)) ImixCn)) 


0 0.99S9g943E+01 

1 0.77879934E+01 

2 D,60fa52795E+01 

3 0.47236280E+01 

4 0.3678748&E+01 

5 0.28&49957E+01 

S 0 .22312472E+01 
7 0 .1737&839E + 01 

3 0.135329a2E+Oi 

3 0.10539396E+01 

) 0.82080042E+00 

. 0 .&3923275E+00 

! 0.49782866E+00 

1 0.38770416E+00 

0.30193993E+00 
0.23514734E+00 
0. 18312991E+00 
0»23514734E+00 
0.30193993E+00 
0.38770416E+00 
0. 497S28S6E+00 
0. 63923275E+00 
0.82080042E+00 
0. 10539396E+01 
0. 13532982E+0i 
0. 1737S339E+01 
0.22312472E+01 
0.23649957E+01 
0.3&73748&E+01 
0.4723G280E+01 
0.60&52795E+01 
0.77a79934E+01 


0.99999998E-02 0 
0 • 1557r3007E-01 0 
0 . 18195892E-01 0 
0. 188946 I3E-01 0 
0. 18393397E~01 0 
0.17190183E-01 0 
0 . 15618986E-01 0 
0 . 13901768E-0i 0 
0 , 121S0012E-Q1 0 
0.10539747E-01 0 
0.90291636E-02 0 
0.76711616E-02 0 
0-&4721378E-02 0 
0.54232136E-02 0 
0.452943S4E~02 0 
0.37626789E-02 0 
0.31135035E-02 0 
0.37626789E--02 0 
0.45294384E-02 0 
0.542a2136E-02 0 
0.64721378E-02 0 
0.76711&16E-02 0 
0.90291G.36E-02 0 
0. 10539747E-01 0 
0. 12i80012E-01 0 
0. 1390176SE-01 0 
0. 156ia98bE--01 0 
0.17190133E-01 0 
0. 18393897E-01 0 
0.18894613E-01 0 
0.13195B92E-01 0 
0 . ir5576007E-01 0 


t s a slowly 
A b s ( X ( n ) ) 


rt S '3 3 ' j *1'* 3 3 i ) E 1 

. 7 7 3 8 0 0 9 2 E T 0 1 

. 6 0 6 5 3 0 6 2 E + 0 1 

.47236657E4-01 

. 38787946 E-^01 

.28650472E+01 

.22313018E+01 

.17377394E+01 

,13533530E+01 

.10539923E4-01' 

.32085001E+00 

.&3927877E-t-00 

.49787074E+00 

.3a774216E+00 

.30197391E+00 

.23517744E-5-00 

. 13315&37E+00 

.2:3517744E + 00 

.30197391E+00 

« 3877 42 1 b E 0 0 

»49787074E+00 

.63927877E+00 

. 82085001E+00 

, 10539923E+01 

..13533530E + 01 

. 17377394E+01 

.22313013E+01 

.28650472E+01 

.3&787946E+01 

.47236657E+01 

.60853062E+01 

.77380092E+01 
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POWER SPECTRAL 


TA 


DENSITY QE 


3, .12 

THE CLUTTER FROM ONE RANGE 


CELL 


( THE PSD VALUES ARE EVALUATED AS FRACTIONS OF THE TOTAL 
POWER RETURNED FROM THAT RANGE CELL ) 


FREQUENCY 

POWER 


CUMULATIVE 

INDEX 

(NORMALISED) 

p in dB' 

POWER <dB) 

( m ) 

C p) 


(ONE-SIDED) 

0 

0«47713768E-i'00 

-■0.32135630E + 01 

-0.18978031E+0i 

1 

0« 197b5860E-^00 

“0 « 70403435E-^01 

-0.39251748E+00 

2 

0,40718611E-01 

-0. 1390207 lE+02 

-0.r3305725E + 00 

3 

0.1 2673 114E-01 

-0. 1397il65E+02 

-0 .&1803705E-01 

4 

0.4476S20&E-02 

-0.234903G3E+02 

-0.35190552E~01 

5 

0.22375344E-02 

-0.2&502304E->-02 

-0.21947678E-01 

6 

0.11187978E-02 

-0,295124a7E+02 

-0. 15340543S-Oi 

7 

0.7243476&E-03 

-0.31400530E-5-02 

-0.11068582E-01 

3 

0.44260218E-03 

-0.33539864E+02 

-0. 346031 12E-02 

9 

0.33673967E-03 

-0. 347264 14E+02 

-0.64767990E-02 

10 

0.23546igiE-03 

-0 .36230792E+02 

- 0 « 5 0 9 0 6 0 5 0 E - 0 2 

11 

0.20il6336E-03 

-0 .36964516E+02 

-0.3906S400E~02 

12 

0. 15578169E--03 

-0.38074833E+02 

-0.29899587E-02 

13 

0. 1459380 9E-03 

-0 .38356827E+02 

-0 . 2 130933 IE -O'" 

14 

0. 12315245E-03 

-0.39095570E+02 

-0. 14066435E-'' 

15 

0, 125121 19E-03 

-0 .39 026695 E +02 

-0.67096431H- 

16 

0.11409923E-03 

-0 .39427177E+02 

0.000 0 00 OOE 

17 

0. 12512716E-03 

- 0 ■ 3 9 0 2 6 4 8 5 E + 0 2 


18 

0. 12315485E-03 

-0.39095490E+02 


19 

0. 1459S399E-03 

-0.33356804E+02 


20 

0. 1557S221E-03 

-0. 3807432 lE+02 


21 

0.2011635iE-03 

-0.3S964503E+02 


oo 

XU uj 

0.23546156E-03 

-0.36280300E+02 


23 

0.3367B749E-03 

-0.34726440E+02 


24 

0.4426021SE-03 

-0.335393&4E+02 


25 

0.7243488SE-03 

~0.31400524E+02 


26 

0. lil87982E-02 

-0.29512434E+02 


■7 7 

0.22375348E-02 

-0,.26502304E + 02 


28 

0.44763234E-02 

-0. 2349030 1E+ 02 


29 

0. 12673126E-01 

-0. 18971165E+02 


30 

0 . 40718656E-01 

-0 . 13902066E + 02 


31 

0. 19765893E+00 

-0..70408363E+01 
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( 
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0 0, 47713766E+00 

1 0. 19765873E+00 

2 0.40718&15E-01 

3 0, 12&7:3118E-01 

4 0. 4476820 lE-02 

5 0.22375358E-02 

6 0. 111S7992E~Q2 

7 0.72434935E-03 

8 0.44260250E-03 

9 0.33679037E-03 

10 0 .23546259E~03 

11 0.20116477E-03 

12 0.15578243E-03 

13 0 . 14598986E-03 

14 0.12315555E-03 

15 0. 12512758E-03 

16 0.11409910E-03 

17 0.12512758E-03 

18 0.12315555E-03 

19 0.14598986E-03 

20 0.15578243E-03 

21 0-20116439E-03 

22 0.23546265E-03 

23 0.33673993E-03 

24 0.44260229E-03 

25 0.72434394E~03 

26 0.11187989E-02 

27 0.22375367E-02 

28 0.44768206E-02 

29 0. 12673 11 3E-01 

30 0.40718615E-01 

31 0. 197&5363E+00 



This program computes DFT directiVu 
COMPLEX SUM(32),C,G 
DIMENSION Y(32) 

S = 0 . 0 

PI=.31415927E1 

0PEN(UNIT = 1 ,NAME='' CELL. DAT'' , IYPE= ■' OLD ' ') 

DO 700 K=l,32 
SUM(K)=(0. ,0, ) 

L=K“1 
REWIND 1 

DO bOO N~i,32 

M=N--1 

I=LAM/32 

J = L.^H-IA32 

B = FLOAT(.J)API/lb. 

G--= ( 0 . , 0 . ) 

I F ( J , E Q » 0 ) G = ( 1 » j 0 . ) 

IE(J.EQ.8) G=(0.,1.) 
lECJ.EQ.lb) Q=(-1.,0.) 

IE(J.EQ.24) G=(0.,-i.) 

IE(CABS(G> ,EQnO. ) Q=CMPLX ( COS ( B ) .-SIN(B) ) 
READ (1,500) C 
iOO E0RMAT(1X,E15.3) 

SUM(K)=SUM(K>+CAG 
>00 CONTINUE 

'00 CONTINUE 

REWIND 1 
CLOSE ( UN IT=1) 

DO 900 K=l,32 
Y ( K ) = < CABS ( SUM ( K ) ) ) kk2 
300 CONTINUE 

DO 1000 K=l,32 
S-~S + Y(K) 

1000 CONTINUE 

DO 2000 K=l,32 
Y<K')=YCK)/S 
2000 CONTINUE 

OPEN ( UN IT = 2,N.AME='' PSDBF.DAT' , TYPE= '' OLD ' ) 

HQ 4000 K=l,32 
WRITE(2„30Q0) K-1,Y(K) 

WRITE (5, 3000) K-1,Y(K) 

3000 FORMAT (IX, I2,1X,E15.3) 

1000 CONTINUE 
REWIND 2 
CLDSE(UNIT=2) 

STOP 

END 

\ kkkkkkkkkkk PROGRAM END kkkkkkkkkkk 



CHAPTER IV 


CONCLUSION 

4.1 SUMK'^ARY OF WORK DONE 

The aim of this work was to develop computer programs 
in FORTRAN IV which could be used to conduct a preliminary 
investigation into some properties of radar clutter, based on 
the statistical evaluation of measured data. The data is to 
comprise of estimates of the samples of complex reflectivity 
of the source, as a function of both discrete range and discrete 
time, when illuminated by a stationary beam. Three properties 
were identified for computer-aided investigation, namely, 
probability distribution, correlation and power spectral density. 
Three distributional models (Rayleigh, '^eibull and Lognormal) 
are commonly suggested for clutter data. It would be of interest 
to investigate the adequacy of these models. The correlation 
between the clutter variations as measured from two range cells 
which are separated spatially, is also of interest. Based on 
the assumption of wide-sense stationary, uncorrelated scattering 
(WSSUS), the scattering function description of a specific range 
cell is the PSD, The estimation of clutter PSD from the data 
is once again of interest. Three programs were developed during 
the course of this work: 

Proqram-I i 

With this program the clutter data can be tested for the 
adequacy of three types of probability distributions, nam.ely the 
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Rayleigh, Weibull and Lognormal, which are commonly suggested 
as valid models. The method employed is a qualitative one. 

The resultant data is required to be plotted. The plot is a 
straight line if the particular model is adequate. The 
program processes clutter data pertaining to one range cell. 

An optional facility for time averaging the data initially is 
provided. 

Proqram-II t 

This program provides two methods for evaluating the 
cross correlation between the clutter returns obtained from 
two range cells: 

a) Direct method 

b) FFT method 

The particular method is selected by the user. As a result of 
program design, in the direct method, the number of data points 
over which correlation is to be evaluated, the start-point in 
the first sequence and the time-delay in the second sequence 
are user-selected variable parameters, and the value of cross- 
covariance (normalised) is computed. In the FFT method on the 
other hand, the value of cross-correlation (normalised) is 
computed over the full length of the input sequences, for all 
values of the delay. In the second method, the sequences are 
assumed to be periodic. The FFT algorithm used is the same 


as in Program-Ill 



72 


Program-Ill ; 

The program evaluates power spectral density (normalised) 
based on the DFT of the clutter data from one range cell. A 
radix 2 in-place decimation-in-time FFT algorithm is used in 
computing the DFT, 

Constructed data was required, to test the programs. A 
program was, therefore, also written to generate ^ayleigh, 

Weibull and Lognormal distributed numbers. Programs I, II and 
III were run with generated test data. They were found to carry 
out their processing operations as designed, 

4o2 SOME SUGGESTED IMPROVEMENTS 

A scheme for preprocessing the raw measured data, i,e. 
one that essentially re-arranges the data cell-wise, needs 
software implementation. It involves transposition of a matrix 
of large size (e.g. 1024x300 size, each element in turn consisting 
of two real numbers denoting a complex sample). (Section 1,2 .3 
refers). 

The test for the adequacy of the distributional models 
is based on the visual assessment of a graphical plot and is 
qualitative in nature. Therefore, the inclusion of a second 
test, which is quantitative, is suggested. One such test is 
the chi-squared test [4] (^ection 2,4 refers). 
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APPENDIX A ; COMPUTER PROGRAMS 

A.l Program-I : An aid to the verification of Hayleigh, 

V/eibull and Lognormal models for the 
clutter data from one range cell. 

A. 2 Program-II : Estimation of correlation between the 

clutter returns from two range cells. 

A. 3 Program-Ill : Estimation of Power Spectral Density of the 

clutter from one range cell. 

A. 4 A program for the generation of Rayleigh, Weibull and 
Lognormal distributed numbers. 

A, 5 List of data files required while using the programs. 
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V A ^ A A A ^ A A A A ^ A A A A A A A A A ^ A A A ^ ^ A A A A ^ A A A A ^ A ^ A A A ^ A A A ^ A A ^ ^ A ?*? A A 

PROGRAMS EGR STATISTICAL EyALUATION OF MEASURED CLUTTER DATA 


M « T ^ L. 1 i o P i' O i O b ;i X S C‘ Sj ^ -_i -U V -3 S W '3 X :f i.% 9 :Ii E L* 5? p *b ^ X X JL A 3 Pi Q U .P It 


P R 0 Q R A M - ' 


PROGRAM - II 


A d e q u a c y o f R A Y L E I G H , U B I B U L L and L 0 Q N 0 R H A L 
Models for the clutter data from one Range 
Cell - A Q u a 1 i t a t i v e A s s e s s m e n t „ 

CORRELATION between the clutter data from 
Two Range Cel 1 s . 


PROGRAM - III 

A A A A A A A A A; A A A. AAA A A A A A A A A A A A A A A A A A; A A A A A A. A A A A. A A A A A A A A A A A A. A A A A. A A A A A A A A A A A A A : 


POWER SPECTRAL DENSITY of the clutter from 
one Range Call based on DEI of the data. 
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k^'kit'h'kh^:khAA^'k:Khk'A'kiii^k'k^\hAJi'kAA'kkA'kx'kiik'kh^iAh^(?i'h^‘.'khhk'k'k:kAk kkkkkkkkkkkkk 


PSOGSAM - r 


THIS PSOQRAM ENABLES A QUALITAIiyE ASSESSMENT OE THE ADEQUACY OF THE 
FOLLOWING DISTRIBUTIONAL MODELS FOR THE CLUTTER DATA 


(A) 

RAYLEIGH 

(B) 

WEI3ULL 

(C> 

LQGNORMA 


ikkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 
PROGRAM MODELS 
COMPLEX UAL 

DIMENSION AZ(I024) .BZ ( 1024 ) ,CZ ( 1024) ,DZ(1024) 

COMMON A2/AREA1/BZ/AREA2/C2/AREA3/DZ 
WRITE (5,2) 

E0RMAT</////////15X, 'ENTER THE NO. OF (COMPLEX) DATA POINTS'// 

• 1 15X,'T0 BE OPERATED UPON FROM FILE CELL. DAT'// 

2 15X,'( NOT TO EXCEED 1024 ) : ',$) 

READ(5,3,ERR=1) NUM 
F0RMAT(I4> 

IF((NUM.GT.1024).0R. (NUM.LT.D) GOTO 1 
JJ=NUM 

0PEN(UNIT=:1 ,NAME='CELL.DAT' ,TYPE='OLD.' ) 

OPEN (UN IT=2, NAMES 'AVERAG.DAT' ,TYPE='OLD' ) 

DO 5 Jsl,NUM 
READ(1,4) UAL 

4 F0RMAT(1X,E15.8) 

WRITE<2,4) UAL 

5 CONTINUE 
REWIND 1 
REWIND 2 
CL0SE(UNIT=1) 

CL0SE(UNIT=2) 

6 WRITE(5,7) 

7 F0RMAT(///36X, 'MENU'/3&X, ' '///> 

WRITE(5,10) 

10 FORMAT(20X, '1 ' ,&X, 'TEST FOR RAYLEIGH/WEIBULL MODEL'// 

1 20X, '2' ,6X, 'TEST FOR LOGNORMAL MODEL'// 

2 20X, '3' ,6X, 'AVERAGE THE DATA OVER GIVEN TIME PERIOD' / 

3 20X,'4' ,&X,'EXIT'/// 

4 26X, 'ENTER A NUMBER : ',$) 

20 READ(5,30,ERR=40) I 

30 FORMAT(Il) 

IF(I.GT.4) GOTO 40 
IF(I.EQ.3) CALL I IMEAV ( NUM » J J ) 

IF(I.EQ.l) CALL WEITST(JJ> 

IF(I.EQ.2) CALL LNTEST(JJ) 

IF(I.EQ.4> STOP 
GOTO & 

40 URITE(5,50) 

50 FORMAK/ZEfaX, ' INCORRECT ENTRY. ENTER AGAIN : ',$) 

GOTO 20 

STOP 'END OF RUN' 

END 


kkk SUBR-l : FINDS LINEAR TIME AVERAGE OF THE ORIGINAL DATA kkk 



SUBROUIINE TIMEAMtN.J) 76 

COHPLEX CC<1024) ,SUM,C • 

URITE<5,2) 

EOSMAT(/i6X,'L = NO. OF CONSECUTIVE DATA VALUES TO BE AVERAGED' 
1 //16X, 'ENTER L 
SEAD<5,5) L 
FORMAK 14) 

IE<L.GT.N) GOTO 90 

OPEN < UN IT=1,NAME=' CELL. DAT ' ,TYPE=' OLD' ) 

J=0 

K=1 

LL=L 

LLL*L 

0 SUM=<0.,0.) 

j=a+i 

DO 30 1=1, LLL 
READ(1,20) C 
10 E0RMAT(1X,E15.8) 

SUM=SUM+C 
)0 CONTINUE 

CC<a>=SU«/ELOATCLLL> 

IF<a.EQ.1024) GOTO 70 

K=K+L 

LL=LL+L 

IF<K.GI.N) GOTO 70 
IFCLL.GT.N) GOTO 60 
GOTO 10 

30 LLL=N+L-LL 

GOTO 10 

70 REWIND 1 . 

CL0SE(UNII=1) 

0PSN(UNIT=2,NAME='AUEEAG.DAT',TYPE='0LD' ) ■ 

URITE(2.20) <CC< I) , I=1,J) 

REWIND 2 
CL0SE(UNIT=2) 

WRITE<5,80) J 

30 F0RMAT<//21X, 'AVERAGED DATA HAS ',14,' COMPLEX SAMPLES'//) 

GOTO 200 

90 URITE(5,100) N 

100 FQRMAT(24X,'L EXCEEDS TOTAL NO. OF DATA <=',14,' )'//) 

GOTO 1 
200 RETURN 

END 


kkk SUBR-2 ; CARRIES OUT STEPS TO CHECK THE ADEQUACY OF THE 
kkk RAYLEIGH AND WEIBULL MODELS.- kkk 

SUBROUTINE WEITSKNO) 

COMPLEX B(1024) 

DIMENSION R1 < 1024 ), El < 1024 ),X( 1024) ,Y<1024) 

COMMON R1/AREA1/F1/AREA2/X/AREA3/Y 
URITE<5,5) 

5 F0RMAT<25X, 'WEIBULL DISTRIBUTION TEST RESULTS') 

WRITE(5.8) 

) F0RMAT<25X,' ') 

OPEN (UN IT= 1 ,NAME=' AVERAG.DAT' ,IYPE='OLD' ) 

READCl.lO) (3(1), 1 = 1, NO) 

.0 FORMAK IX, E15. a; 

DO 00 1 = 1, NO 
Rl( I)=CAE3(B( I) ) 

\{S OiTKiT T'*.i I IT7 
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WRITE<5,30)AME,VAR 

F0aMAT<23X, "MEAN=' ,E17 . 3/29X , ' VAR IANCE= ' ,E17 . 8/ ) 

CALL CUM(NQ,NOE) 

DO 40 J-=1,MQE 

Y< J)=AL08(AL0Q( (J) 5 ) > 

CONTINUE 

URITE<5,45> 

FORMAT (3 IX, 'X' , ISX , " Y V31X , 16X , ' - ' ) 

URITE<5^50) (X(L) ,Y(L) ,L=1,NQE) 

. FQRMAT<23X.2E17.8) 

WRITE<5,52) 

I F0RMAT(/31X, 'PLOT <X,Y) POINTS') 

WRITE<5,55) 

5 F0RMAT<27X, ' IF POINTS LIE ON A 3TR LINE') 

WRITE(5,57) 

? • F0RMAT(26X, 'DISTRIBUTION IS WEI3ULL W ITH ' /29X . ' ETA=SLOPE ' ) 
CALL EIT<NO,NOE,SL,YIN) 

A=YIN/SL 

SIG=EXP(-A) 

WRITE(5iS0)SIQ 

0 E0RMAT(29X, 'SIGMA=' ,E17.8/) 

URITE(5,70) 

0 FORMAKSSX, 'BEST STR LINE FIT IS :') 

URITE(5,80> SL.YIN 

0 FORMAKSgX, 'SL6 pE=' ,E17.8/29X,'Y INTERCEPT*' , E17 . 8/ > 

REWIND 1 
CL0SE<UNIT=1) 

0PEN<UN1T=2,NAMH = 'UBTEST.DAT' ,TYPE*'OLD' ) 

WRITE(2,90)<X<I) ,Y( I), I=1,N0E) 

0 F0RMAT(1X,2E15.8) 

REWIND 2 
CL0SE(UNIT=2) 

RETURN 

END 


*** SUBR-3 : ARRANGES MAGNITUDE VALUES OF THE DATA IN ASCENDING 

ORDER. AAA 

SUBROUTINE ASC(N) 

DIMENSION R(1024) 

COMMON R 
M = N-1 

DO 20 1<=I,M 

DO 10 L=K,M 

IF(R(K) .Le.R(L+l)) GOTO 10 
T=R(K) 

R(K)=R<L+1 ) 

R<L+1)=I 
0 CONTINUE 

5 CONTINUE 

IF(R(N> .NE.0.0) GOTO 40 
WRITE<5,30) 

) F0RMAT(25X, 'ALL DATA ARE ZERO ' //25X, 'STOPP INQ' ) 

STOP 

) RETURN 

END 


AAA SUBR-4 : EVALUATES MEAN AND VARIANCE AAA 
SUBROUTINE MV(N,AM.V) 
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p=6 

Q=0 

C-=N 

DO 10 1=1, N 
P=P+R( I> 

CONTINUE 

AH=P/C 

DO 20 J=1.N 

Q=Q+<R( a)-AM>A*2 

CONTINUE 

V=Q/C 

RETURN 

END 


Ah SUBR-5 : FINDS THE VALUE OF THE CUMULATIVE DISTRIBUTION 
AA FUNCTION AT EACH DATA POINT. 

• SUBROUTINE CUM<N,NS> 

DIMENSION R(1024) ,F < 1024 ), EX< 1024 ) 

COMMON R/AREA1/F/AREA2/EX 

IS=0 

NS=N 

DO 10 K=1,N 
IF(K.EQ.N) GOTO 20 

IE« (R<K) .NE.R(K+1)).AND.(R(K) .NE.0.0)) GOTO 5 

NS=NS-1 

GOTO 10 

IS=IS+1 

E< IS >=FLOAT(K) /FLOAT <N) 

EX( IS)=ALOG<R(K) > 

CONTINUE 
F(NS)=. 999999 
EX(NS>=ALQG(R<N) ) 

RETURN 

END 


AAA 8UBR-6 : FITS BEST STRAIGHT LINE THROUGH (X,Y) POINTS 
AAA IN MAIN PROGRAM 

SUBROUTINE F IT < N ,NS , S, YT) 

DIMENSION A(1024) ,B(1024) 

COMMON /AEEA2/A/AREA3/B 

Q=NS 

ASUM=0 

BSUM=0 

AASUM=0 

ABSUM=0 

DO 100 1=1, NS 

ASUM=ASUM+A( I) 

BSUM=BSUM+B< I) . 

AASUM=AASUM+A( I)AA2 
ABSUM=ABSUM+A< I)AB< I) 

• CONTINUE 

D=GAAASUM~ASUMAA2 

S=<GAABSUM-ASUMABSUM)/D 

YT = < AASUMABSUM“ASUMAABSUM>/D 

RETURN 

END 


AAA 

AAA 


AAA 

AAA 


ril tn 


AAA LC6N0SMAL MODEL. AAA •? 

SUBROUTINE LNTESKNO) ' 

COMPLEX B(1024) 

DIMENSION SS(1024) . EF ( 1024 ) .X ( 1024) , Y ( 1024) 

COMMON RR/AREA1/EE/ARSA2/X/AREA3/Y 
URITE<5,5) 

EQRMAX(23X, 'LOGNORMAL DISTRIBUTION TEST RESULTS' /23X, 

1 / 

0PEN(UNIT=1 .NAME= 'AVERAQ.DAT' ,TYPE='OLD' ) 
READ(1,10)(B(I),I=1,N0> 

E0RMAT(1X,E15.8) 

DO 20 1=1, NO 
RR( I)=CABS{B(I)) 

» CONTINUE 

CALL ASC(NO) 

CALL CUM(NQ,NOE) ' 

DO 30 1=1, NO 
Y( I)=RR< I) 

RR( I)=ALOG(RR( I) ) 

> CONTINUE 

CALL MV<NO,AME.yAR) 

URITE(5,40) AME,VAR 

) EQRMAT(28X, 'MEAN=' ,El&.a/28X,'UARIANCE=' .E?w6.3/28X, 

1 'OF LN OF DATA'/) 

DO 50 1=1. NO , 

RR(I)=Y(I) 

) CONTINUE 

S=SQRT(UAR) 

llIPTTF^®^ 

5 F0RMAT<35X,'X<,16X,'Y'/35X,'-',16X.'-') 

DO SO J=1,NQE 
X< J)=(X( j)-AME)/S 
Y( J)=EE( J> 

A=ABS(X(a) ) 

CALL ERECA.E) 

IF(X<J)) 60,70,70 . 

0 X(J)=.5-E 

GOTO 80 
X( J)=.5+E 

«RITE<5,85) J,X(J).Y(J) 

F0RMAT(18X, I4,4X,2E17.8) 

CONTINUE 
«RITE(5,g5) 

5 F0RMAT(/32X. 'PLOT (X,Y) PO INIS ' /27X , ' IF POINTS LIE ON A STR LIN 

1 /28X, 'DISTRIBUTION IS LOGNORMAL'/) 

CALL EIT(NO,NOE,SLO,YI) 

WRITE(5,100) SLO.YI 

00 E0RMAI(30X,'BEST STR LINE FIT IS:'/ 

1 30X,'SL0PE=' ,E16.8/27X,'Y INTERCEPT=' ,E1 6. 8 ) 

PLOT (X< I) ,Y( I) , 1=1 ,NOE) 

PLOT Y=SLOAX+YI 
REWIND 1 
CL0SE(UNIT=1) 

0PEN<UNIT = 2.NAMS='LNTEST.DAT'.TYPE='0LD' ) 

WRITE (2. 200) (X( I), Yd) ,1 = 1. NOE) 

00 EaSMAT(lX,2S15.8) 

REWIND 2 
CLOSE ( UN :T=2) 

RETURN 



«: SUBR-8 rCOMPUTES THE VALUE OE THE ERROR FUNCTION kAk 
SUBROUTINE ERE(H,Q) 

DIMENSION U(16),W(16) 

F(X)=EXP(-X*X) 

DATA U<1) ,U(2) ,U(3) ,U<4>/.04a3077, .1444720,.2392874, .3313686. 

DATA U<5>,U(6),U<7) ,U( 8)/. 4213513, .5068999, .5377153, .6630443/ 

DATA U<9) ,U( 10) ,U( 11) ,U<12)/.7321821, .7944833, .8493676, .8963212/ 
DATA U(13),U(14),U(15),U( 16) /. 9349061 , .9647623, .9856115, .9972639/ 
DATA U<1),W<2),W<3) ,W < 4 )/ .09654009 ,. 09563872 , .09384440 . .09117388/ 
DATA W<5) ,W(6) ,U<7),W(8)/.08765209..08331192,.07819390, .07234579/ 
DATA W(9) ,W(10) ,W(ll)/.06582222,.05868409,.0509g806/ 

DATA U(12) ,W(13) ,W( 14)/. 04283590, .03427386, .02539207/ 

DATA W<15) ,W(16)/.01627439,. 00701861/ 

Q=0 • 

IF(H.EQ.O. ) GOTO 20 
A=H/(2.*SQRT(2. ) ) 

DO 10 I--l,16 
X*A*U< I) 

. X1*A+X 
X2=A'-X 

Q=Q+U( I)*(F<X1)+F(X2) ) 

CONTINUE 
Q=Q*AA. 5641896 
RETURN 

END \ kkkkkkkkkkkk END OF PROGRAM - I ************! 
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PROGSAH - II 


THIS PROGRAM EVALUATES THE 
CLUTTER RETURNED FROM TWO 
OF TWO METHODS PROVIDED : 

CORRELATION 
RANGE CELLS. 

< SPATIAL ) BETWEEN 
THE USSR CAN SELECT 

THE 

EITHER 

<A) 

DIRECT 

METHOD 


<B) 

FFI 

METHOD 



kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk 
PROGRAM CORREL 
COMPLEX SEQ(1024) 

COMMON /AREAl/SEQ 
WRITE(5,100) 

FORMAT (//////////////////////) 

URITE<5,200) 

FORMAT(///6X, 'CHOOSE A METHOD ;-'//) 

WRITE<5,300) 

FORi1AT(23X,'l DIRECT'// 

23X,'2 EFT'///) 

URITE(5,400) 

FORMATtSX, 'ENTER 1 OR 2'// 

&X,'( Any a1;hsr character EXITS )'// 

6X, 'ENTRY 

READ(5,500,ERR=600) M 
EORMAT< 12) 

IE(<M.NE.1).AND.<M.NE.2>) GOTO SOO 
WRITE(5.100) 

CALL NUMDAT(NUM) 

WRITE(5,100) 

lE(M.EQ.l) CALL CORDIR<NUM) 

IE(M.EQ.2) CALL CORFFKNUM) 

GOTO 50 

STOP 'END OF RUN' 

END 


kkk SUBR-l: COMPUTES CORRELATION COEFFICIENT ( COVARIANCE, NO.. »ISED 
kkk BY THE DIRECT METHOD , USING THE VALUES OF 'NOS*, * I3IART 

kkk AND *ITAU* SUPPLIED TO IT. 

SUBROUTINE CORDIRtNO) 

COMPLEX DAT,DATK<.1024) ,DATKK(1024) ,UK,UKK,SUM,RO 
SUM=(0. ,0. ) 

WRITE(5,5) 

S FQRMAT(////////////////31X, 'C ORRELATIO N'/31X. 

I / //) 

CALL TIMES ( NO , NOS , ISTART , ITAU ) 


OPEN (UNIT=1, NAME* 'ONESET. DAT' ,TYPE=' OLD') 
IF(ISTART.EQ.l) GOTO 30 
M=ISTART-1 
DO 20 1=1. M 
READd.lO) DAT 
10 FORMATdX.ElS.S) 

20 CONTINUE 

30 CONTINUE 

DO 40 3=1. NOS 


50 

00 

00 

00 

1 

.00 

1 

2 

iOO 


SCO 
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lO 


30 


60 

70 


80 


90 


100 


200 


READd.IO) DAXK-d) 

SUM=SUiirl!ATK{ J) 

CONTINUE 

UKR=SEAL< SUM) /FLOAT < NOS) 

UK I=AIMAG( SUM) /FLOAT (NOS) 
yK=CM?LX(UKS,UKI) 

DO 50 1=1, NOS 
DATK( I)=DATK(.I)-IJK 
CONTINUE 
REWIND 1 
CLQSE(UNIT=1) 

SUM=(0. ,0. ) 

0PEN(UNIT=2.NAME='TU0SET.DAT' ,TVPE='OLD' ) 

IF( ( ISTART .EQ . 1) . AND . ( ITAU.Sa.O) ) GOTO 70 

M=ISTART+ITAU-1 

DO 60 1=1, M 

READ(2,10) DAT 

CONTINUE 

CONTINUE 

DO SO J=1,N0S 

READ<2,10) DATKK(J) 

SUM=SUM+DATKK(3) 

CONTINUE 

UKKR=aEALC SUM) /FLOAT (NOS) 

UKKI=AIMAQ< SUM) /FLOAT (NOS) 
UKK=CM?LX(UKKR,UKKI) 

DO 90 1=1, NOS 

DATKK( I)=DATKK( I)-UKK 

CONTINUE 

REWIND 2 

CL0SE(UNIT=2) 

SUM=(0„ ,0. ) 

SUMKD=0 


SUMKKD=0 
DO 100 J=1,N0S 

SUMKD=SUMKD+(CABS(DATK( J))) A^2 
SUMKKD=3UMKKD+(CABS(DATKK(3) ) )*jk2 
CONTINUE 

DENOM=SQST<SUMKD)*SQ&T(SUMKKD) 

DO 200 1=1, NOS 
A=REAL(DAIK( I) ) 

AA=REAL(DATKK( I) ) 

B=AIMAG(DATK( I) ) 

BB=-AIMAQ(DATKK( I) ) 

C=AAAA-BABB 


D=AABB+AAAB 

SUM=SUM+CMPLK(C,D) 

CONTINUE 

PI=4.AATAN(1. ) 

CF=130./PI 

EOS=REAL(SUM)/DENQM 

RQI=AIMAG(SUM)/DENOM 


ROMAQ=CABS ( SUM >/DENOH 
IF(ROR.SQ.O.O) R0PHAS=90. 


IE( (SOS.EQ.0.0) .AND. ( (SDI/RGMAG) .LT.0.0) ) SOPHAS=-ROPHAS 
;E(EQMA3.EQ.O.O) R0PHAE=0. ! Ip. this case ROPHAS is not defipe 
I F ( S ■ 1 M S .. 0 . 0 ;■ R S ? H A S = A T A M ( A & S ( R 0 1 / S C R ) ) A C E 
[F( (ROR.NE.0.0) .A---D. ( (ROI/RQR) .LT.O.C) ) RDPHAS=-R0PHA3 
>iRITE(5,300> RQR.ROI.RCMAQ.SOPHAS 
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26X, 
2£X, 'r 


RETURN 

END 


r = '.E15,S' + )'// 

3bX, ' I r I •= %EI3.8// 

26X,'< r ( in deg. ) = SS15.8////) 


^kk 

SUBR-2: ASKS 

kkk 

(a) 

kkk 

(b) 

kkk 

(c) 

kkk 

IT F 


SUBROUTINE I 

WRITE(5.10) 

F0RMAT(//14X 


iTART* ). kAk 

(b) = ‘HIAU*. AAh 


1 OVER WHICH'/14X, 'THE CORRELATION COEFFICIENT IS TO BE COMPUTED 

2 > : '//14X, "CORRELATION LENGTH* ',$) 

READ(5,20,E8S=25) NOSAMP 
FORMAT ( 14 :> 

IF(NOSAMP.GE.l ) GOTO- 40 
WRITE(5,30) 

FORMAT(//28X, " INCORRECT ENTRY :ENTES AGAIN") 

GOTO 5 

IF(NOSAMP.LE.N) GOTO 60 
URITE(5,50) N 

FaRMAT(//15X, "CORRELATION LENGTH EXCEEDS NO. OF INPUT DATA (=", 
1 14,')') 

GOTO 5 
WRITE<5,70) 

E0RMAT<//18X, "ENTER TAU (TIME-DELAY) AS THE NO. OF INTERVENING" 
1 /18X, 'DATA POINTS : ' //18X. ' TAU=' , $ ) 

READ<5,20,E8R=75) HTAU 
IF<MTAU.6E,0) GOTO 90 
WRITE<5.80) 

FORMAT (//30X, ' INCORRECT ENTRY (TAU) 

GOTO 60 

IF( (MTAU+NOSAMP) .LE.N) GOTO 200 
WRITE(5,100) N 

E0RMAT(//26X, " (TAU+CORRELATION LENGTH) EXCEEDS" 

1 /26X," INPUT DATA LENGTH ( = ",14,")') 

GOTO 60 
WRITE(5.300) 

F0EMAT(//17X, 'THE CORRELATION CALCULATION CAN BEGIN ANYWHERE' 

X /17X,'IN THE INPUT SEQUENCE. ENTER ;') 

WRITE (5, 350) 

FQRMAT(/17X, 'START AT DATA POINT NQ.(e.g. 1, 5, etc. )=',$) 
READ(5,20.ERR=375) MSTART 
IF(MSTART.GE.l) GOTO 500 
WRITE (5. 400) 

F0RMAT(//17X, " INCORRECT ENTRY (STARTING DATA POINT NO.). ENTER 
GOTO 325 

IF( (MSTART+MXAU+NOSAMP-l ) .LE.N) GOTO 300 
WRITE (5. 600) 

FORMAK //12X, ' (STARTING ?0 INT ) +TAU+ < CORRELAT ION LENGTH)-!" ) 
WRITE (5. 700 M 

FORMAT- i2X. 'EXCEEDS INPUT DATA LENGTH (=',14,'), RE-ENTER 

■j U ... L* ^ 

WRITE -j . '300 ) rt’-ITAxT , MTAU , NOSAMP 


800 



5 F0SHAH/////////////28X, 'STASTING DATA POINT NO , = ' , I4/28X . 

1 'TAU=' . I4/28X, 'CORRELATION LENGTH® 14//// ) 

RETURN 

END 


kkk SUBR-3: ASKS THE USER TO ENTER THE NO. OF ( COMPLEX ) DATA AAA 
kkk POINTS STORED IN THE INPUT DATA FILES 'ONESET.DAT' AAA 

AAA AND 'TWOSET.DAT'. AAA 

SUBROUTINE NUMDAT(N) 

write: 5, 2) 

E0RMAT(/18X, 'ENTER THE NO. OF (COMPLEX) DATA AS ONE OF THESE 
l:'//18X,'2 4 8 IS 32 64 123 256 512 1024'// 

2 26X,'DATA LENGTH TO BE PROCESSED =',$) 

READ<5.3,ERR=5> N 
FORMAK 14) 

NC = 1 

DO 4 NCC=1,10 
NC=NCA2 

IFCN.Sa.NC) GOTO 7 

CONTINUE 

WRITE(5,6) 

F0RMAT</26X, 'NO. INCORRECT: ENTER AGAIN') 

GOTO 1 

WRITE(5,8) N 

F0EMAT(/27X, 'DATA LENGTH® ' , 14, IX . 'BE ING PROCESSEL '//) 

RETURN 

END 


AAA SUBR-4: COMPUTES DFT OF A COMPLEX SEQUENCE USING A RADIX 2 / 

AAA IN-PLACE DECIMATION-IN-TIME EFT ALGORITHM. AA 

SUBROUTINE FFT(N) 

COMPLEX A< 1024) ,T,U,W 
COMMON /AREAl/A 
C=.31415927E1 
DO 5 M=l,10 
IF(M.EQ.IO) GOTO 10 
N1=2AAM 

lE(N-Nl) 5,10,5 
CONTINUE 
0 IM=M 

NN=N-1 
3 = 1 

DO 50 1=1. NN 
IFd.GE.J) GOTO 20 
T=A<J) 

A( J)=A( I) 

A < I ) =T 
0. K=N/2 

0 IF(K.GE.J) GOTO 40 

3 = J-K 
K=K/2 
GOTO 30 
0 J= I + K 

0 CONTINUE 

DO 80 L=l.If1 

U=(l . ,0. ) 

II=2aAL 

T ■? r. T T / •*> 
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IE( I2>QT«2) W=CMPLX(CQS<B) ,-3IN<B) ) 
DO 70 aj=i.:2 

DO 60 K1<=JJ,N,II 

KM=KK+I2 

I=U*A<HH) 

A(MM)=A(KK)-T 

A(KK)=A(KK)+1 

CONTINUE 

U=UAU 

CONTINUE 

CUNTINUE 

RETURN 

END 


AA 

AA‘ 

AA 


SUBR-5: ASSUMES THE TWO INPUT DATA SEQUENCES TO BE PERIODIC 


AND CJIMPUTHS 

FOR ALL TIME 


THE VALUE OF CURSELATION ( NORMALISED ) 
SHIFTS BY THE FFT METHOD. 


SUBROUTINE COSFFTCN) 

COMPLEX B(1024) ,C(1024).D<1024) 
COMMON /AREAl/D 
C0N=45./ATAN(1. ) 

SUMBB=0 
SUMCC=0 
CALL HEADFT 


0PEN(UNIT=1,NAME='0NESET.DAT%TYPE:='QLD' ) 
READ(1,10) (D(I) ,I=1,N) 

F 0 RMAThx,E 15 . 8 ) 

CALL FFT(N) 


DO 20 1=1, N 
B<I)=CQNJG(D(I)) 


SUMBB=SUMBB+(CABS(B< I) > ) Aa2 

CONTINUE 

REWIND 1 

CL0SE(UNIT=1) 

OPEN (UN IT=2,NAME= 'TWOSET.DAT' ,TYPE='OLD' ) 
READ(2,10) (D( I),I = 1,N) ■ 

CALL FFT(N) 

DO 30 2=1, N 
C< J)=D( J> 

SUMCC=SUMCC+(CABS<C<J) ) )AA2 

CONTINUE 

REWIND 2 

CL0SE(UNIT=2) 

DEN=SQRT<SUMBB)ASQRT(SUMCC) 

DO 40 1=1, N 
D< I)=(B< I)AC( I) )/DEN 
CONTINUE 
CALL FFT(N) 

0PSN<UNIT = 3,NAME = 'C0RFFT.DAT' ,TYPE="OLD' ) 

DO 60 1=1, N 

L=I-1 


RE=REAL<D( I) ) 
AI=-AIMAQ(D( I) ) 


ABSL=CAB3(D( I) ) 

IF(RE.EQ.O.O) PH=90, 

IF ( ( RE . EG . 0 . 0 ) . AND . ( v A 1/ A3SL ) . IT . 0 .0 ) ) PH= -PH 

IF(ABSL. gQ.O.O) P;-:=0. ! In this case PH is not defined. 


; '-. A 


t. 


i-i CO ^ 4fj 
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E0RMAT<4X, I4,2X,4(2X,E15.8> ) 

WRITE<3,55) L,RE,AI,ABSL,PH 

EORMAKIX, I4,4( 1X.E15.8> > 

CONTINUE 

REWIND 3 

CLOSE(UNIT=3) 

RETURN 

END 


:* SUBR-6 : DISPLAYS HEADINGS. IT FORMS PART OF THE EFT METHOD. kkA 
SUBROUTINE HHADFT 
WRITE(5,1) 

EORMAK 15X, 'CELL-TO-CELL CROSS CORRELATION'/ 

/ / //lOX 

'r = (COMPLEX-VALUED) CORRELATION COEFE IC lENT '//3X , 
'Time-del av' ,55X, 'Phase (r ) '/3X, 

' Index ',&X,'Re(r)',12X.'Im(r)',12X.'Abs<r)',9X,' indeg.'/3X, 

/ S12X,' ' ,12X,' S9X,' ' ) 

RETURN 

END [ kkkkkkkkkkkk END OF PROGRAM - II kkkkkkkkkkkk\ 
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PSQQSAfI - III 


THIS PROGRAM EVALUATES THE POWER SPECTRAL DENSITY OF THE CLUTTER 
FROM ONE RANGE CELL, GIVEN THE TIME-DOMAIN SEQUENCE CONSISTING 
OF THE ESTIMATES OF THE COMPLEX REFLECTIVITY OF THAT RANGE CELL. 


AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAC 
PROGRAM PSD 


00 


}0 


00 


too 


1 


500 


1 

100 

'00 


'50 

iOO 

00 


COMPLEX B(1024) 

DIMENSION BB(1024) ,DD(513) .CUM (513) 

COMMON /AREAl/B 

SUM=0.0 

SUM2=0.0 

CALL NUMDAT(NUM) 

CALL HEADER 

0PEN<UNIT=1 ,NAME="CELL.DAT' , TYPE= ' OLD " > 
READ(1,100> (B( I) .I = 1,NUM) 
F0RMAT(1X,E15.8) 

REWIND 1 
CL0SE(UNIT=1) 

CALL FFT(NUM) 

DO 200 3=1, NUM 
BB(3)=(CABS(B<a)))AA2 
SUM=SUM+BB( J) 

CONTINUE 

NNN=NUM/2+l 


DD<1)=2.ABB<1) 
DO 300 M=2,NNN 
IARQ=NUM-M+2 


DD(M)=BB(M)+BB( lARQ) 

CONTINUE 

SUM2=SUM+BB<1)+BB(NNN> 

CUM(1)=DD(1) 

DO 400 N=2,NNN 
aARQ=N-l 

CUM(N)=DD(N>+CUM( JARG) 

CONTINUE 

0PEN<UNIT=2,NAME='PSD.DAT' ,TYPE='OLD'> 

DO 700 L=1,NNN 
LLL=L-1 

WRITE <2, 500) LLL, (BB(L)/SUM) , <10.AAL0G10(BB(L)/SUM) > , 
<10.AAL0G10(CUM(L)/SUM2) > 

EORMATdX, I4,3(1X,E15.8) ) 

WRITE (5, 600) LLL, (BB<L)/SUM) , (10.AALOG10(BB(L)/SUM) ) , 
(10.AAL0Q10(CUM(L)/SUM2) ) 
E0RMAT(16X,I4.,3(2X.E15.8) ) 

CONTINUE 

MMM=NNN+1 


DO 900 aK=MMM,NUM 
LLL=JK-1 

WRITE( 2,750) LLL , < BB < JK ) /SUM ) , ( 10. AALQQIO < BB ( JK ) /SUM ) ) 
FOKMAKIX, I4,2(lX.E15.a) ) 

WRITE (3,800) LLL, <3B(3K)/SUM) ,(10. AALOG10(BS(aK)/SUM) ) 
?0 R M AT ( 1 6X , 14 . 2 ( 2X , E i 5 . S > ) 

CGNTINUE 

rs i.j t m n -5 


n r) o o o rj n o 



STOP 'HMD CE RUN 
END 
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hiiA SUBR-l: GENERATES HEADINGS ON SCREEN 
SUBROUTINE HEADER 
WRITE(5,10) 

0 EORMAK//////////////////////) 

WRITE<5,20) 

Q. FORNAK 12X POWER SPECTRAL DENSITY OF THE CLUTTER FROM' 
1 ' ONE RANGE CELL' 

■ 2 

3 / 

WRITE(5,30) 

>0 FORMAT( 12X. ' < THE PSD VALUES ARE EVALUATED AS FRACTIONS 

1 'OF THE TOTAL' 

2 /12X, 'POWER RETURNED FROM THAT RANGE CELL )'//) 
WRITE <5. 40) 

10 FORMAT (12X, 'FREQUENCY POWER 

1 'CUMULATIVE' 

2 /12X,' INDEX (NORMALISED) p in dB 

. 3 'POWER (dB)') 

WRITE(5,50) 

30 F0RMAT(/12X.' (m) (p) 

1 '(ONE-SIDED)' 

2 /12X,' 

3 /) 

RETURN 

END 


9 


/ 


SUBR-2: ASKS THE USER FOR THE NO. OF < COMPLEX ) DATA POINTS TO M 
khk BE OPERATED UPON FROM FILE CELL.DAT. ki 

SUBROUTINE NUMDAT(N) 

WRITE(5,1) 

EORMATdSX, 'ENTER THE NO, OF (COMPLEX) DATA AS ONE OF THESE :'// i 
1 13X.'2' 4 8 16 32 64 128 256 512 1024') i 

». READ(5,3) N i 

I FORMAT (14) 

NC=1 . 

DO 4 NCC=1,10 

IE(N.eLnC) goto 6 
\ CONTINUE 

WR ITE ( 5 5 ) 

S FaRMAT(/28X,'N0. INCORRECT: ENTER AGAIN') 

GOTO 2 ; 

i WRITE(5,7) N 

' F0RMAT(/27X, 'DATA LENGTH^' .14, IX, 'BEING PROCESSED'///) 

RETURN 

END i 

kkk SUBR-3: COMPUTES THE DFT OF A COMPLEX SEQUENCE USING A RADIX 2 kk-, 
kkk IN-PLACE DEC It; AT ION- IN-TIME FFT ALGORITHM kk) 

SUBROUTINE FFT(N) 

COMPLEX A(1024) ,T.U,W i 

COMMON ^ARHAl/A 
C= .31415-32721 
DO 5 M=1.13 
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N1=2AAM 

lE(N-Nl) 5,10,5 

CONTINUE 

IM=M 

NN=N-1 

J=1 

DO 50 1=1, NN 
lEd.GE.J) GOTO 20 
T=A(J) 

A< J>=A< I) 

A( I)=T 
• K=N/2 

» IF<K.GE.a) GOTO 40 

J=J-K 
K = K/2 
GOTO 30 
> J = a+K 

) CONTINUE 

DO 80 L=1,IM 
U=(1.,0.) 

II=2AAL 
12=11/2 
B=C/ELOAI( 12) 

IE(I2.EQ.l) W=<-1.,0.) 

IE(I2.EQ,2) W=(0.,-1.) 

IE( I2.QT.2) W=CHPLX(COS(B) ,-SIN(B) > 

DO 70 aa=i,i2 
RRspEALCJ) 

GQ=AIMAi3<U) 

DO 60 KK=JJ,N. II 
MH=KK+I2 
R=REAL( A(MM) ) 
a=AIHAG(A<MM) ) 

TR=R*RR-GAGG 

TQ=RAQQ+RR*8 

I=CMPLX(TR,TQ) 

A(MM)=A(KK)-T 
A(KK)=A<KK)+T 
0 CONTINUE 

R=REAL<W> 

Q=AIMAG<W) 

XR=R^RR-GAGQ 

TG=R*GG+RRAQ 

U«UAW 

U=CMPLX(TR,TQ> 

0 CONTINUE 

0 CONTINUE 

RETURN 

END [kitMAMk END OF PSCORAH - III kkkkkkkkkkl 



4 


90 


i'ki(Jtkki^:k^i^-k^kkkkkkkk^k^-kAAkkkkkk^:kkk^kA:kkkkhk^^k^ik^:k*ii^kAAkAkkA^kkkkC 

C 


PROGRAM FOR TEST DATA SIMULATION 


C 


THIS PROGRAM GENERATES : 

<A) WEIBULL 

CS> LOGNORMAL 

DISTRIBUTED REAL NUMBERS ARRANGED IN ASCENDING ORDER 
OF MAGNITUDE. 


C 

C 

C 

c 

n 

L: 

C 

p 

w 


kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkC 
PROGRAM RANDOM 
WRITE(5,2) 


FORMAT (//////////////////////) 

WRITE(5,3) 

FORMATOIX, ^NUMBER GENERATOR V31X, . 

1 

WRITE (5.4) 

E0RMAT<////22X, ' 1 WEIBULL DISTRIBUTED NUMBERS V/22X, 

1 '2 LOGNORMAL DISTRIBUTED NUMBERS V///22X, 

2 • 'ENTER 1 OR 2 < ANY OTHER NO. EXITS ) 
READ(5,5,ERR=G) I 

FORMAKIl) 

IFd.EQ.l) CALL WBL 

IF<I.EQ.2) CALL LOQNML 

IF( < I.NH.l) .AND.(I.NE.2) ) GOTO 6 

GOTO 1 

& STOP 'END OF RUN' 

END 


*** SUBR-l: ASKS THE USER TO ENTER THE VALUES OF THE DISTRIBUTION 
kkk PARAMETERS 'ETA' AND 'SIGMA' DESIRED, AND GENERATES 

kkk REAL NUMBERS, IN ASCENDING ORDER OF MAGNITUDE, WITH A 

kkk WEIBULL/RAYLEIQH DISTRIBUTION. 

SUBROUTINE WBL 
DIMENSION F<1024) 

CALL NUMIiAT(N) 

WRITE<5,1) 

FORMAK lOX, 'ENTER WEIBULL PARAMETERS ETA, SIGMA.'/ 

1 10X,'( ETA = 2.0 FOR RAYLEIGH )'// 

2 lOX. 'ETA=' ,$) 

READ (5,2) ETA 
F0RMAI(E15.a) 

WRITE<5,3) 

FORMAK lOX, 'BIQMA=' , t) 

READ(5,2> SIGMA 
FORMAK IX, E15. 3; 

FGRMAT(27X, I4.3X,E15.3) 

M=N~1 

00 10 1=1. M 

Sil) = : IC-AI V I /FLOAT •; ;• 

CONTINUE 


kkk 

kkk 

kkk 

kkk 
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T = 1./(1."F( I) ) 

U=ALOG(T) 

V=:ALCQ(U>/ETA 

W=V+ALOG<SIQHA) 

X=EX?(y) 

«RITE(5,5) I,X 

WRITE(1,4) X 

CONTINUE 

REUIND 1 

CLOSE(UNIT=l) 

RETURN 

END 


*** SUBR-2: ASKS HOW MANY NUMBERS ARE TO BE GENERATED FROM THE USER AA? 
SUBROUTINE NUHBAT(NUM) 

WRITE(5,1) 

FORMAT C 16X, "HOW MANY NUMBERS ARE TO BE GENERATED ?'// 

1 IGX, 'ENTER IT < SHOULD NOT EXCEED 1024 ) : ',$) 

READ<5,3,ERR=4> NUM 
F0RMAT<I4) 

IF( <NUM.LT. 1) .OR. <NUM.GT.1024> ) GOTO 4 

GOTO 6 

WRITE<5,5) 

F0RMA1</26X, 'NO. INCORRECT; ENTER AGAIN') 

GOTO 2 

WRITE<5,7> NUM 

F0RMAT</24X, 'DATA LENGTH*' , 14, IX, 'BEING PROCESSED'///////) 

RETURN 

END 


AAA SUBR-3: COMPUTES THE VALUE OF THE ERROR FUNCTION AAA 
SUBROUTINE ERE(H.Q) 

DOUBLE PRECISION X1,X2,2 
DIMENSION U<16),M(1&) 

E<Z)=DEXP<-ZAZ> 

DATA U< 1) ,U<2) ,U< 3) ,U<4)/. 0483077. .1444720, .2392874, .3318&8G/ 

DATA U(5) ,U<&) ,U( 7) ,U<8)/. 4213513, .5068999, .5877158 ,. 6630443/ 

DATA U(9) ,U<10),U<11),U(12)/.7321821, .7g44838,.8493676. .8963212/ 
DATA U(13) ,U(14) ,U<15) ,U< 16 )/ .9349061 , .9647623 , .98561 15 , .9972639/ 
DATA U( 1) ,W<2) ,W<3) ,W<4)/.0g6540O9, .09563872, .09384440, .09117388/ 
DATA U<5) ,U<6) ,W(7) ,U(8)/. 08765209, .08331192, .07819390, .07234579/ 
DATA W ( 9) ,W< 10), well)/. 06582222, .05868409,. 05099806/ 

DATA W(12) ,W(13) ,W( 14)/. 04283590, .03427386, .02539207/ 

DATA W(15) ,W( 16)/. 01627439,. 00701861/ 

Q=0 

IF(H.EQ.O. ) GOTO 30 
A=ABS(H)/(2.ASQRT<2.)) 

DO 10 I«l,16 
X*AAU( I) 

X1=A+X 

X2=A-X 

Q = Q<-W< I)A<F(X1)+F(X2) ) 

CONTINUE 

a=aAA 

I?<H) '20,30,30 
Q--' 

'SSTL'RN 
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SU£t]k"-4; L''n;.ER FOS^THE PES*lISSIELE ERROR ( ABSOLUTE VALUE ) 

^ _ LOGNOEKAL -D ISTR IBUTED NUMBERS AND 
IriE VrtuUES QE iHE DISTRIBUTION PARAMETERS 'MU' AND 
'SIGMA' IiESIRED. 

SUBROUTINE PERM<ERE ,UU,SIQMA) 

HRITE(5,22) 


FORMAT < IbX, 'ENTER THE VALUE OE ABSOLUTE ERROR PERMISSIBLE'/ 
15X,'IN THE COMPUTATION AS . 

READ(5,24,ERE=20> ERR 

FORMAT(F10.7) 

«RITE(5,26) 

E0RMAT(/////21X, 'ENTER LOGNORMAL DISTRIBUTION PARAMETERS'/ 
21X,'MU AND SIGMA AS — - •'// 

1X,'MU=' ) 

READ<5,28,ERR=25) UU 
FORMAT (F7. 3) 

WRITE<5,30) 

F0RMAT(/1X, 'SIGMA=' ) 

READ(5,32,ERR=2S) SIGMA 
FORMAT (F7. 3) 

WRITE(5,34) 

FORMAT < ///////////35X. ' GENERATED DATA' /35X , 

'//) 

RETURN 

END 


:k 

kA A 
AAA 


SUBE-5: GENERATES REAL NMBERS IN ASCENDING ORDER OF MAGNITUDE 

WITH A LOGNORMAL DISTRIBUTION. AAA 

SUBROUTINE LOGNML 
CALL NIJMDAT(NO) 

P«NO 

CALL 'PERM<PE,U,SIG) 

PE=PEA1. 7724538 

OPEN(UNIT = l,NAME='LOG.DAT',IYPE = 'OLri' ) 

X=1.7724538/P-. 8862269 
CALL ERF(X,T) 

DO 90 1=1, NO 
F=<FLOAT( I>/P)A1. 7724538 
IFd.EG. (N0/2> ) GOTO 30 
IFd.EQ.NO) GOTO 40 
Y=.8862269+T“F 
IE<ABS(Y) .LE.PE) GOTO 50 
C=.l ■ 

YO=y 

<o=x 

ro=T 

<=XO“CAyO/ABS<YO) 

:all erf(x,t> 

'».8862269+T-F 

:F<ABS(Y> .LE.PE) GOTO 50 

:F<ABS(Y) .LT.ABS(YC) ) GOTO 10 

:=xo 

=T0 
= . 1 A C 

FCC.LT. ( 1 .E-iO) ) GOTO 50 
OTC 20 
= 0 . 
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1 = 0 . 

GOTO 50 

40 X=4. 534099 

50 Z=EXP(XaS1G+U) 

60 WRITE (5, 70) I,Z 

70 F0RMAT(27X, I4,3X,E15.8) 

WRITE (1,80) Z 
80 FGRMAT(lX,E15-8) 

90 CONTINUE 

REWIND 1 
CL0SE(UNIT=1> 

RETURN 

END END OF TEST DATA SIMULATION PROGRAM 



94 


A.5 List of data files required while using the programs. 

All input, output and intermediate data files used in any 
of the programs are required to be created before run-time. The 
maximum length of data which the programs are designed to handle 
is a one— dimensional array of 1024 (complex) sample values. 
Program-wise list of data files is given below: 

Proqram-I 


Program Name : PROGRAM MODELS 
Input data file 


File name : CELL, DAT 

Nature of data: Complex reflectivity samples of one range cell 

Type of data : complex 

^ata format : IX, E15.8/lX,El5.8 

(Each complex sample is stored as two successive real numbers, 
being its in-phase and quadrature components, to be immediately 
followed by the next complex sample in the sequence). 


Data entry mode z by the user 
Intermediate storage file; 

File-name ; AVERAG.DAT 

Nature of data: Time-averaged values of input data 

Tyne of data : complex 

i^ata format : 1X,E15.8/1^,E15.8 

Data entry mode : by program 

(Time-averaging is an optional facility provided in the program) 



Output data files: 

a) File-name ; WBTEST.DAT 
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Nature of data: Set of (x,y) points to be plotted as a result 

of the ^eibull/Rayleigh test program segment 
processing, (Each record in the output 
consists of two real numbers corresponding 
to the values of X and Y respectively). 

Data format : IX, 2E15.8 

Data entry mode: by program 

b) File -name : LNTEST.DAT 

Nature of data : ^et of (x,y) points to be plotted as a result 

of the Lognormal test program segment 
processing. (Each record in the output 
consists of two real numbers corresponding 
to the values of X and Y respectively). 

Data format : IX, 2E15.8 

Data entry mode: by program 


Proqram-II 


Program name : PROGRAM CORREL 

Input data files: 

a) File-name : ONESET.DAT 

Nature of data' : Complex reflectivity samples of the 

first range cell. 

Type of data : Complex 

Data format :, IX, E15.8/1X,E15.8 

Data entry mode : by the user. 
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b) File-natoe : TWOSET.DAT 

Nature of data : Complex reflectivity samples of the 

second range cell. 

Type of data : Complex 

Data format : IX, E15.8/1X, E15.8 

Data entry mode : by the user 

Intermediate storage files : none 
Output data file (for FFT method only): 

File-name : CORFFT.DAT 

Nature of data : Set of normalised correlation values, 

corresponding to all values of time-delay 
by the FFT method. The exact output consists; 
of a set of 5 real numbers per record, 
corresponding to 

m; Re(r);Im(r)} Abs(r)j Ph(r) (deg) 
where 

m is the time delay index 

r is the (complex-valued) correlation 
(normalised) and 

Re(.), Im(.), Abs(.) and Ph(.) 

denote respectively its real part, 

imaginary part, magnitude and phase. 

Data format : IX, 14, 4(1X, E15.8) 

Data entry mode : by program. 

Note: Output in case of the direct method is a single value of 

normalised cross covariance, which is displayed on the VDU, 
and no output data file is required. 
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Program-Ill 

Program name : 

Input data files: 
File-name 
Nature of data 

Type of data 
Data format 
Data entry mode 
Intermediate storage 
Output data file: 
File-name : 

Nature of data : 


Data format : 
Data entry mode : 


PROGRAM PSD 

: CELL, DAT 

: Complex reflectivity samples from one range 
cell. 

: Complex 

: IX, E15.8/1X, E15.8 
: by the user 
file : none 

PSD. DAT 

Normalised PSD values corresponding to all 
values of the frequency index (k). The exact 
output consists of four real numbers per 
record corresponding to k; PSD (normalised) (p)i, 
p(in dB) j cumulative value of PSD (one-sided) 

(in dB) . 

IX, 14, 3(1X, E15.8) 
by program 



Program for the generation of Rayleigh, Weibull and 
Lognormal distributed numbers: 


^8 


Program name 
Input data files 
Intermediate storage 
Output data files: 

a) File-name 

Nature of d ata 


I^ata format 
i^ata entry mode 

b) *‘ile-name 

Nature of data 


^ata format 
^ata entry mode 


PROGRAM RANDOM . 
none 

files : none 

: WBL.DAT 

: Set of real numbers arranged in ascending I 
order of magnitude, which satisfy a 
Rayleigh or Weibull distribution, depending 
upon the parameter values entered by the 
user. 

: IX, E15.8 

: by program 

; LOG. DAT 

: Set of real numbers arranged in ascending 
order of magnitude, which satisfy a 
Lognormal distribution, in accordance 
with the parameter values entered by the 
users 

: IX, E15.8 

; by program 
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Note 1: Output data files (a) and (b) are required for storing, 
respectively, the generated (a) Weibull/Rayleigh and 
(b) Lognormal distributed numbers. Only one of these 
need be created, if data is to be generated pertaining 
to the concerned distribution only. 

2. Two further processing steps are required before the data 
as generated above can be used as test data for Program-I, 
viz. (a) jumbling the generated numbers in pseudo-random 
fashion, (b) creating complex numbers, whose magnitude 
values are equal to the numbers in (a) (e.g. by assigning 
a phase value to each in pseudo-random fashion). 
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